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SUMMARY 


An analysis is presented for Earth-orbiting spacecraft missions designed 
to measure stratospheric constituents using solar occultation. This technique 
uses the Sun as an energy source and depends on the attenuation of solar radia- 
tion as the Sun's image rises or sets through the atmosphere as viewed from a 
spacecraft. The orbital mechanics of solar occultation are examined, and some 
results are summarized for a simulated 1-year mission using an orbit with an 
inclination of 70° and an altitude of 600 km. Of particular concern is the 
spatial and temporal atmospheric sampling relative to the Earth's surface. The 
nominal orbit has been chosen so that some spatial coverage can be obtained at 
all latitudes. There are 1 to 15 opportunities per year, depending on latitude, 
to make a series of consecutive measurements in the vicinity of any particular 
latitude. Sampling capabilities of solar occultation missions are then examined 
parametrically for a range of orbit conditions. 

A simulated 1-year data base is presented for the nominal mission. The 
data consists of 9480 "measurements" of a single quantity associated with each 
measurement opportunity. These data are produced by a model whose underlying 
functional variability is not revealed. The measurements are grouped into lat- 
itude bands, and the observed variability within latitude bands is assumed to 
be both systematic and random. The statistical properties of the data set are 
examined to determine the extent to which global behavior can be extracted from 
solar occultation mission data. A standard curve-fitting technique is applied 
to these data, by which parameters are determined for a simple periodic mathe- 
matical model with variable mean, amplitude, and phase. 


INTRODUCTION 

Global monitoring of the stratosphere is important in attempting to under- 
stand and protect the Earth's environment. Stratospheric components of partic- 
ular interest are aerosols, ozone, and many trace gases which participate in 
the complex photochemical processes of ozone production and depletion. These 
components, through their ability to scatter and absorb solar radiation, play a 
critical role in the Earth's radiation balance and hence in maintaining climate 
and weather within acceptable limits. 

The use of Earth- or biting spacecraft to provide reliable long-term global 
monitoring is part of a total commitment to providing the necessary data base 
for present and future environmental studies. One technique for measuring 
stratospheric constituents using a spaceborne sensing system is to measure the 
attenuation of solar radiation as the Sun's image rises or sets through the 
atmosphere. There are two monitoring opportunities per orbit, one sunrise and 
one sunset, except for those times when the spacecraft is in continuous line-of- 
sight contact with the Sun. These so-called solar occultation measurements are 
potentially applicable to a variety of stratospheric monitoring activities. In 



one recently proposed flight experiment, the objective is to learn more about 
the distribution of stratospheric aerosols (ref. 1); in another proposal, the 
objective is to study halogens (such as HCl ) in the stratosphere (ref. 2). 

In the former case, the proposed flight experiment attempts to expand a body 
of knowledge that is already quite large, to the extent that many facets of 
aerosol production, distribution, transport, and depletion are understood 
reasonably well. In the latter case, the importance of halogens in strato- 
spheric ozone chemistry has only recently been recognized on theoretical 
grounds, and is unverified by direct measurements. Hence, the proposed flight 
experiment is only a first step toward understanding concentration levels and 
the coarsest features of the underlying distributions. A shared characteristic 
of these missions, typical of cases where solar occultation measurements are 
expected to be useful, is the need for global long-term coverage. This paper 
demonstrates how this basic requirement can be met and provides a preliminary 
outline of one way in which the resulting data might be used to determine the 
behavior of stratospheric constituents over spatial and temporal scales of 
interest. 

The concept of ‘’global" coverage is a flexible one which may be inter- 
preted in different ways for meeting different objectives. For the present 
purposes, global coverage is interpreted to mean that all latitudes are covered 
by a single satellite on one or more occasions during a year. This requirement 
leads to orbits of high, but not necessarily polar, inclination, as shown when 
the orbit geometry is examined in detail. It is unfortunate that the global 
coverage requirement cannot be met with Sun-synchronous orbits so that solar 
occultation measurements would be compatible with many other types of Earth- 
monitoring activities, but this specialized geometry gives only restricted lat- 
itudinal coverage over a range which depends on the initial timing of the orbit 
with respect to the Sun. Such a case has been examined in a previous study 
(ref. 1), where the objective was to show how a two-satellite mission could 
effectively combine the coverage of a Sun-synchronous orbit with that achieved 
at a lower inclination. 

An inherent limitation of the solar occultation technique is that the 
local time of day (the "clock” time) at the locations of the measurements is 
necessarily restricted to dawn or dusk. Thus, the technique is not capable of 
monitoring diurnal variations. These must either be assumed to be insignifi- 
cant or be accounted for with mathematical modeling and/or data from some other 
source. 

The overall possibilities for and limitations on spatial and temporal cov- 
erage with the solar occultation technique must be assessed on the basis of the 
total amount of data obtained in a specific period of time. To achieve some 
understanding of how such data would be accumulated, a twofold approach has 
been used. First, a nominal 1-year mission has been defined and its particular 
output examined in detail. This mission has been selected to achieve some 
global coverage in the sense described previously. Then, from the insight 
provided by a specific case, some important aspects of the nominal mission are 
generalized and expanded parametrically to cover a wider range of orbit options. 
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After the orbital analyses are complete, the nominal mission is recalled 
to generate a set of simulated measurements for statistical examination. For 
the purposes of this paper, it has been assumed that adequacy of spacecraft 
coverage can be investigated in terms of the extent to which spatially and tem- 
porally varying distributions can be extracted from a simulated set of single- 
parameter values. These data are produced by a single-valued function of 
longitude, latitude, and time which associates one scalar quantity (the 
"measurement") with each occultation opportunity. The actual process of 
recording a complete sunrise or sunset is much more complicated; it consists 
of a series of individual measurements obtained as the Sun rises or sets, each 
of which must then be interpreted in terms of variations in the measured quan- 
tity. One such interpretation, but certainly not the only one, would be to 
assume that the measurements constitute a vertical profile through the atmo- 
sphere. This involves data reduction algorithms which are of no concern for 
the present analysis. Further treatment of the data, such as integration over 
altitude, could yield a single parameter associated with each sunrise or sun- 
set. Integrating a vertical profile could result in a total vertical burden, 
although other single-parameter results might also be of interest. However the 
data analysis is accomplished, another basic assumption of this paper is that 
the details are of no concern for assessing the adequacy of global coverage of 
the occultation events. 

A third assumption about the simulated data is that for the proposed 
investigation of coverage capabilities with a particular spacecraft-based 
measurement technique, the specific origin and form of the single-parameter 
model is not important. It can produce any set of values whose spatial and 
temporal variability includes and exceeds that which can be extracted from the 
resulting data base. In this paper the model is thought of as a "black box" 
incorporated into a computer program for orbit propagation. At each sunrise 
or sunset, the appropriate space and time coordinates are input to the model, 
which returns a single value. It might be desirable to have this scheme pro- 
duce output which an experimenter can relate to his expectations and past expe- 
rience, but this is not necessary. The details of the model used in this paper 
are not discussed at all, but are revealed only through the measurements simu- 
lated in the nominal mission. This approach is an accurate representation of 
the real environment wherein underlying distributions of observed quantities 
are revealed piecemeal through the measurements which are available. Any other 
approach would place unwarranted emphasis on the suitability of the model for 
predicting the actual behavior of a physical quantity - a goal which is totally 
beyond the scope of this study. 

In the statistical analysis of the data set, two different procedures are 
followed. First, some descriptive statistical parameters are calculated, mean 
values, for example, and confidence limits are discussed for this type of 
information obtained with no knowledge of or assumptions about the functional 
form of the variability in the data. Second, some assumptions are made about 
the data in the form of a simple periodic model, and three parameters - mean, 
amplitude, and phase - are fitted to the data. Confidence in this process is 
then related to the amount of the total observed variability which can be 
represented by the model, without regard to its physical basis. 
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SYMBOLS 


a semimajor axis, km 

A amplitude, one of three parameters in a model for q, units of q 

B coefficient matrix in the model for q 

B,b,c,d angles, defined in figure 10, needed to calculate Xqi given i, 

h, and 5 q, deg 

C matrix of parameters , €2^ and C3 to be estimated by the least 

squares technique in the model for q 

D central angle between a spacecraft and the Earth’s horizon, deg 

f true anomaly, deg 

h altitude, km 

unit heading vector (see fig. 2) 

inclination, deg; also, a summation index 

longitude, deg 

matrix of observations of q 

number of values in a sample for statistical analysis 

unit vector perpendicular to a spacecraft orbit (see fig. 2) 

point on the surface of a spherical Earth at which solar occultation 
measurement is made 

single parameter assumed to have been extracted from a series of mea- 
surements taken at sunrise or sunset, expressed in arbitrary units 

sample mean calculated from a set of q’s, units of q 

Q. mean, one of three parameters in a model for q, units of q 

r spacecraft position vector from the center of the Earth, km 

rqi position vector of tangent point from center of Earth 

t0 radius of Earth, 6378.145 km 

R ratio which determines how well variability in data has been 

represented 
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RA right ascension, deg 

s^ sample variance, units of 

S vector in direction of X 0 defined in equation ( 1 ) 

SC spacecraft 

t time, days, or dimensionless ratio normalized to 1 year 

Y velocity, km/sec 

Sc 0 unit vector to the Sun 

Y an arbitrary observed quantity 

a, 3 pointing angles to the Sun referenced to a spacecraft coordinate 

system (II,N,f) , deg 

Y phase factor, one of three parameters in a model for q, deg 

5 declination, deg 

e residual error in a model for q, units of q 

Ti angle defined as 90° + 60 (see eq. (5)) , deg 

X latitude, deg 

H actual mean (unknown) of a population whose mean is to be estimated 

from observations made on the population, units of q 

T direction of the vernal equinox vector, for locating the x-axis of 

an Earth-equatorial (right-ascension— declination) coordinate system 

0 ) argument of perigee, deg 

right ascension of ascending node, deg 

Subscripts (not defined in main symbol list) : 

rel relative to the horizon at the time of a sunrise or sunset 

R residuals, as when calculating the sum of squares of the residuals 

T tangent point at the surface of a fictitious spherical Earth or on 

a concentric spherical shell around the Earth 

TG tangent point referenced to the surface of an ellipsoidal Earth 

0 referring to the Sun 
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Notation: 


dimensioned vector 

unit vector; or to denote estimates of the value of scalar quantities 

% 

averaged quantity obtained from a set of observations 


THE ORBIT DYNAMICS OF SOLAR OCCULTATION MISSIONS 

The well-known precession of a spacecraft orbit plane, as predicted by 
first-order perturbation theory (ref. 3) , is of dominant interest for the 
analysis of solar occultation missions. Both the argument of perigee O) and 
right ascension of the ascending node ^ are functions of time. This changing 
Sun— orbit-plane geometry determines the limits on and rates of spatial coverage. 
The geometry for characterizing a sunrise or sunset viewed from Earth orbit 
is defined in figure 1. All the vectors are expressed in a right-ascension- 
declination (RA- 6 ) system, which is an Earth-equatorial system with the x-axis 
in the direction of the vernal equinox vector (T) . The unit vector xq points 
from the spacecraft to the center of the Sun. The orbit plane is located in 
space by its inclination to the equator i and by the angular elements ca) 
and The spacecraft is located within an orbit by its true anomaly f, the 

angular distance between perigee and spacecraft. A sunrise or sunset occurs, 
by definition, at that instant at which the projection of X 0 is tangent to 
the surface of a spherical Earth at Prp; the position vector to Pip is trjx, 
with longitude and latitude coordinates L.p and Xip. (This differs from the 
usual astronomical definition of sunrise or sunset as viewed from the ground, 
but it is an appropriate choice for this study. It can be modified if necessary 
to account for the oblateness of the Earth and other contributions to a varia- 
tion in the actual location of the horizon as viewed from the spacecraft.) 

The first step in locating the tangent points is outlined in appendix A. 
(This appendix has previously been published as appendix C of ref. 3; it is 
included here for completeness.) A derivation is given for finding the orbit 
true anomalies at which occultations occur. Given such a point on the orbit, 
its radius vector r (in the RA -6 system) can be found (see, for example, 
ref. 3) . Then, from simple geometrical considerations, rrp is the vector sum 
of r and a vector which is in the direction of X 0 and whose magnitude 

is 
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(1) 


The derivation in appendix A assumes fixed orbit elements. Since both U) and 
9, are changing, some iteration is required. Fortunately, the changes in oj 
and are only fractions of a degree per orbit. So, the procedure which has 
been followed is to calculate the true anomalies at which occultations occur 
for an orbit whose angular elements oa and 9, are fixed at the values 
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corresponding to the time of perigee. Then, the times from perigee to the 
points of occultation are used to update w and and the calculations are 
repeated. This one-step iteration gives occultation points to a sufficiently 
high degree of accuracy for the intended purpose. 

Note that the maximum latitude of the tangent point can exceed the maximum 
subsatellite latitude and that the location of the tangent point depends simul- 
taneously on the orbit inclination and altitude and on solar position. It is 
useful to think of the tangent latitude as depending predominantly on the orbit 
parameters and the subsolar declination and the tangent longitude as depending 
largely on the Earth's rotation. Consider the angular rates involved: rela- 
tive to an inertial coordinate system, the Sun appears to process at about 
1 deg/day in the positive direction; an orbit of less than 90° inclination 
processes in the negative direction at a typical rate of a few degrees per day, 
while the Earth rotates in the positive direction at 361 deg/day. Thus, it is 
expected that the tangent latitude would undergo a relatively slow variation 
as the Sun's declination changes and the orbital plane precesses, while the 
Earth's rotation would produce a much more rapid variation in the tangent 
longitude. 

Turning now to the definition of the vector X0 relative to a spacecraft, 
figure 2 establishes a coordinate system at the spacecraft for defining point- 
ing angles to the Sun; a is a "pitch" angle (down to the horizon for a sun- 
rise or sunset) and 3 is a "yaw" angle which locates the Sun relative to the 
plane of the spacecraft motion. For a circular orbit and a spherical Earth, 
ot has a constant negative value at the horizon for each orbital altitude. 


Analysis of a Solar Occultation Mission 

With the discussion and definitions of the previous section, a hypotheti- 
cal solar occultation mission can be formulated and a set of sunrise and sunset 
data examined. An orbit has been selected with an inclination of 70° and an 
altitude of 600 km. Some relevant data are given in table 1 . This choice of 
orbit altitude and inclination was based on the fact that for moderate alti- 
tudes of a few hundred kiloneters (to ensure lifetimes of several years) , 70° 
is an inclination which allows coverage at each pole sometime during a year; 
thus the orbit provides global coverage as defined earlier. (This matter is 
treated parametrically in the following section.) The orbit was started near 
the summer solstice of 1981 in an orientation (dependent on the time of day 
of injection into orbit) which does not produce occultations for the first 
few days; the initial geometry is shown in figure 3. The orbit so defined is 
henceforth referred to as the nominal orbit. A mission starting with this ini- 
tial orbit configuration and lasting for 1 year is referred to as the nominal 
mission. The launch year is arbitrary and the results could be reproduced for 
any other year desired. 

Figure 4 shows the variation of several important mission parameters as a 
function of time for 5440 revolutions (1 year) of the nominal orbit. Through- 
out this figure (and similar figures throughout the paper) the small crosses 
refer to sunsets, while the dots refer to sunrises. Figure 4(a) shows the 
variation in the tangent latitude, as previously defined in figure 1 , as a 
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function of time from launch. The cycles of coverage are driven by the orbit 
precession relative to the Sun and are typical of nonsynchronous orbits. The 
period of the cycles is about 100 days. This can be verified by noting that 
the precession rate of this orbit is about -2.5 deg/day relative to inertial 
space and, consequently, about -3.5 deg/day relative to the Sun. Note that 
curves in figure 4 are not continuous; they are made up of discrete events 
occurring at the rate of about 30 per day. The slow change in the northern and 
southern latitude limits of the coverage is due to the yearly cycle in subsolar 
declination and some additional orbit constraints. At least the initial dis- 
tribution of sunrises and sunsets can be understood qualitatively by consider- 
ing the rotation of the orbit plane from its initial orientation shown in fig- 
ure 3. For the nominal orbit there are about 9500 measurement opportunities 
in a year, equally divided between sunrises and sunsets. The exact number is 
dependent on the precise timing of the launch, which can be altered to shift 
the coverage cycles within the limits imposed by solar position. It is pos- 
sible to select a launch time (an initial orbit plane orientation) which will 
not permit tangent latitude coverage all the way to the poles, so that some 
planning is necessary in that regard. However, because launch time variations 
do not have a major effect on the overall performance of a 1-year mission, they 
are not investigated parametrically in this study. 

Figure 4(b) shows, for the nominal mission, the variation in pointing 
angle g required to track the sunrises and sunsets throughout the year. It 
is not unexpected that a full ±180° of scanning is required to accomplish this 
task. Note the symmetry evident between the sunrise and sunset curves, which 
may be of considerable significance in setting up a strategy for pointing the 
spacecraft sensors in preparation for detecting a sunrise based on the angular 
position of the previous sunset. 

The tangent longitude coverage is examined separately for sunrises and 
sunsets in figures 4(c) and (d) , respectively. The overall impression is one 
of reasonably complete longitude-latitude coverage over a year, but the way in 
which spatial coverage is built up as a function of time is obscured in these 
two parts of figure 4. This matter is dealt with in more detail in subsequent 
figures . 

Figures 4(a) to (d) treat all occultation events as equivalent and each 
tangent point as a single set of two coordinates. In reality, the changing 
conditions under which these events occur produce variable sunrise and sunset 
rates which affect the amount of time available for measurements. When the 
Sun is in the orbital plane, the apparent vertical motion of the Sun's image 
through the atmosphere is maximum, while at the other extreme, when the space- 
craft is moving parallel to the terminator with the Sun off to one side, the 
sunrise or sunset rate approaches zero. There are potential advantages and 
drawbacks at each extreme which must be taken into account during mission plan- 
ning. One way to keep track of these changes is to calculate the vertical 
motion of the Sun's image relative to the horizon: v^ei = ^SC * ^T *^he 

instant when the projection of xq touches the horizon. When the Sun is 
behind or ahead of the spacecraft, this value is maximum and the length of time 
required for the Sun's image to pass through the atmosphere is short enough to 
consider Vj-ei constant during the measurement opportunity. At the other 
extreme, the motion of the Sun's image in the vicinity of the horizon is more 
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complicated and must be examined in detail. The cyclic variation between fast 
and slow sunrises and sunsets is directly related to the yaw angle 3 , previ- 
ously identified in figure 4(b) . Thus, values of 3 near 0*^ or ±180® place 
the Sun ahead of or behind the spacecraft and result in fast occultations, 
while values near ±90° place the Sun to one side of the spacecraft and result 
in slow occultations. In figure 4(e), Vj-g]^ is plotted for the nominal mis- 
sion. The right ordinate gives the time required for the Sun's image to pass 
through 100 km of atmosphere at the constant rate As noted previously, 

this rate is a good approximation for the larger values of becomes 

less meaningful as the magnitude of ^rel ^PPi^o^ches zero. 

A consequence of the variable measurement time is that the tangent point 
as previously defined moves during the course of a sunrise or sunset. The 
resulting smearing of the surface resolution is of concern for interpreting 
measurements. The correspondence between the time available for measurements 
and the change in surface coordinates is not a simple one, but the limiting 
cases can be described. As the spacecraft reaches the extremes of its subsat- 
ellite latitude, the change in subsatellite latitude with time is near zero; 
the rate of change in tangent latitude is also generally lowest if the occulta- 
tions occur during these portions of the orbit. If the occultations occur when 
the satellite is near the equator, the rate of change of the subsatellite lati- 
tude is maximum, amounting to about (360^/5805 sec) sin 70^ « 0.06 deg/sec 
for the nominal orbit. The tangent latitude changes at about the same rate. 

For the tangent longitude, the situation is reversed, so that the maximum rates 
of change occur at the highest subsatellite latitudes. 

To illustrate the expected changes in surface conditions during a sunrise 
or sunset, some representative measurement situations have been selected from 
the nominal mission. It must be emphasized that selecting specific orbit situ- 
ations for close inspection, as opposed to considering overall mission perfor- 
mance must be viewed with caution. The computer model used to generate occul- 

tation data (relative to a spherical surface) is not intended to actually 
predict the precise occurrence of a particular event many months in advance. 
Even if a spacecraft were maintained in its orbit so that its actual precession 
rate exactly matched the approximate values used in this study, there would 
still be some minor discrepancies in timing of the events because of the dif- 
ference between the actual Earth horizon and the spherical model of the Earth's 
surface. These discrepancies would be most noticeable during the slow sunrises 
or sunsets. However, situations like the ones described in the following para- 
graph are certain to occur during the course of any similar mission. 

A limiting case for a slow sunset and sunrise is examined in detail in 
figure 5. The upper plot shows the altitude of the center of the Sun's image 
above the horizon as a function of time (refraction is ignored) ; the two middle 

plots show the pointing angles a and B required to track the Sun; and the 

lower plot shows the latitude and longitude of the tangent point. For this 
figure, the computer model was temporarily modified to represent the Earth's 
surface as an ellipsoid (see ref. 4) to improve the estimate of where the hori- 
zon actually is. This event is taken from near the beginning of the mission, 
orbit 107 on day 7. The center of the Sun just reaches the horizon and then 
starts back up again, giving more than 12 min for observation of the atmosphere 
below 100 km. On preceding orbits there are other slow passages of the Sun 
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into and through the atmosphere. Table 2 shows the minimum altitude reached 
by the Sun on several consecutive orbits and the length of time during which 
the Sun's center is between the horizon and 100 km. These representative val- 
ues illustrate the possibilities for measurements near the beginning and end of 
periods of occultation. Figure 4(e) shows that there are eight occasions dur- 
ing the nominal mission during which the very slow occultations occur, at which 
time the situations described in figure 5 prevail. Some of the viewing condi- 
tions are summarized in table 3 for the slow sunset and sunrise (orbit 107); 
and they are compared with values for a fast sunset or sunrise taken during the 
period of maximum relative vertical velocity of the Sun's image with respect 
to the horizon (orbit 415). The frequency of viewing conditions between the 
extremes of Vj-gj^ is given in table 4 by actual count and percentage for the 
entire nominal mission. As seen in table 4, nearly all of the occultation 
events occur in less than 100 sec. 

To complete specific analysis of the nominal mission, there remains only 
the additional investigation of tangent longitude-latitude coverage and how 
it is obtained as a function of time. In figures 6(a), (b) , and (c) , the tan- 
gent longitudes and latitudes of sunrises and sunsets are plotted for the first 
7, 30, and 90 days after the start of occultations, respectively. In the first 
of these plots (7 days) the points associated with the first 15 orbits on which 
occultations occur are numbered. It is clear from figure 6 that, typically, 
the 360° range of longitudes is scanned fairly uniformly with a resolution of 
roughly 25° in a time during which the tangent latitude changes relatively 
little. A final view of longitude-latitude coverage is presented in figure 7. 
Here, the Earth's surface has been divided into longitude-latitude boxes of 
approximately equal area, corresponding to a 5° x 5° box at the equator, such 
that the longitude width is an integer factor of 360°. If the latitude band 
width is fixed, the farther a latitude band is from the equator, the greater 
must be the longitude span of each box in that band to give the same area. 

The longitude divisions as a function of latitude are made such that the number 
of boxes within each 5° latitude band results in the minimum departure of the 
area of each box from the reference area at the equator (5° x 5° « 3 x lO^ km^) . 
A summary of this global division is given in table 5 for one hemisphere. It 
shows the latitude band, the area included within each band, the number of 
boxes in the band, the longitude span of each box, the area of each box, the 
cumulative surface area from equator to pole, the percentage of the hemi- 
sphere's area in each band, and the cumulative percentage of hemisphere area 
from equator to pole. Figure 7 shows that the global distribution of tangent 
points is not uniform. There is sparse or no coverage of some boxes in the 
tropics, with heavier coverage toward the poles. (Some parametric examination 
of this point is presented in the next section.) The lack of uniform temporal 
coverage, clearly evident in figure 4(a), presents a potential problem in data 
analysis which is nowhere evident in figure 7. While the ability to obtain 
point measurements in space and time is potentially valuable, it is of more 
interest to be able to separate time and space effects and to determine whether 
the data can be grouped in some way to detect, for example, persistent or per- 
manent spatial differences or trends in some measured quantity over the globe. 
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Parametric Investigation of Tangent Latitude Coverage Available 
From Solar Occultation Missions 

From the spatial coverage capabilities of a specific hypothetical solar 
occultation mission, it is of interest to generalize some of the observed 
results. It is of particular value to understand the factors which influence 
tangent latitude coverage, as this parameter is the most important one for 
determining how much of the globe is covered. For any orbit, the available 
tangent latitudes are constrained to lie within an envelope determined by 
orbit parameters and solar declination. So, first of all, the boundaries of 
the envelope must be determined. Figure 8(a) illustrates, for an orbit of 
inclination i and altitude h, limiting conditions under which a sunrise or 
sunset can be observed, given a solar declination 6. The spacecraft is at its 
maximum subsatellite latitude, and the plane of the figure contains 5^0 and 
N, as defined in figure 2. Generally, orbits would be oriented so that occul- 
tations occur regularly as the spacecraft passes into and out of the Earth's 
shadow. However, for some conditions, as illustrated in figure 8(a), it is 
possible for a spacecraft to be in continuous sunlight. The limiting case^ 
occurs when 


( 2 ) 


angle between spacecraft and tangent point at 


(3) 


the orbit, a function only of the circular orbit 
altitude h. For the conditions shown, sunrises and sunsets occur just at the 
northern extreme of the terminator latitude. (There is an analogous geometry 
for the southern extreme, when 6© is negative.) Such occultations occur, in 
the limit, at Xt = 90° ~ 1<S0|* Note that this is a limit imposed by the lati- 
tude extremes of the solar terminator. Additional constraints may be imposed 
by the ability of a spacecraft to view occultations at any particular time; 
that is, the spacecraft must be able to view the terminator from the dark side, 
looking perpendicular to the terminator from an angular distance D. If h 
or i is increased or if 6© is decreased, in figure 8(a), no occultations 
could be observed for this orientation of the orbit plane. 


D + i + |6©1 =90° 


where D is the Earth central 
the horizon: 


D = cos 


-1 


r© 


r© + h 


The angle D is a constant of 


^If i > 90° (as for Sun-synchronous orbits), use 180° - i in place 
of i for these calculations. 
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Figure 8(b) illustrates a condition under which no periods of continuous 
sunlight ever occur and, as a consequence, the tangent latitude at which an 
occultation event is viewed can never be exactly at the extremes of the termi- 
nator latitude. This occurs when 


D + i + 230.44 < 90° (4) 


where 23®. 44 is the maximum magnitude for the solar declination. The maximum 
tangent latitude, north or south, at which a sunrise or sunset can be observed 
approaches, but may not equal or exceed D + i. 

A third limiting case is illustrated in figure 8(c). Here, for 
D + i = 90°, periods of continuous sunlight can occur regardless of the 
solar declination if the orbit is properly oriented. 

To summarize, 

1. If 230.44 + D + i < 90°, sunrises and sunsets are always observed. 

2. If 1^1 + D + i = 90°, total sunlight is possible for some orienta- 
tions of the orbit plane during part of the year. 

3. If D + i = 90°, total sunlight is possible for some orientations of 
the orbit plane regardless of solar declination, that is, at any time during 
the year. 

The application of these limiting values for solar occultations to a range of 
orbit inclinations and altitudes is summarized in figure 9. Combinations of 
orbit inclination and altitude which fall above a curve of constant solar 
declination permit periods of continuous sunlight for some orientation of the 
orbit plane, while those below the curve do not. Above the 1^1 = 0 curve, 
periods of continuous sunlight may be observed at any time during the year with 
an appropriate orientation of the orbit plane. (The nominal orbit, indicated 
in fig. 9, is such an orbit.) Below the |^1 = 23°. 44 curve, periods of con- 
tinuous sunlight may not be observed regardless of the time of year or orbit 
plane orientation. 

By extension, the data in figure 9 also indicate the maximum tangent lati- 
tudes observable from a particular orbit. To see how this works, consider the 
variation of tangent latitude as a function of orbit parameters. The geometry 
is illustrated in figure 10. The orbit inclination i and solar declination 
^0 are given. The desired quantity is the tangent latitude Xp. Note that 
the angle D is the central angle defined in equation (3) and is constant for 
a given orbit altitude. To solve for the unknown angles in the spherical 
geometry, first let 


T1 = 90° + 


(5) 
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Then, 



Note that when [tan D/tan d[ >1, |sin b| >1 and no occultations can occur. 
These equations hold only for i = 90^. The angle c defines the nodal orien- 
tation of the orbit relative to the solar terminator. For the purposes of 
plotting data from equation ( 6 ) , it is more convenient to use nodal orienta- 
tions relative to the subsolar meridian - RA 0 . Varying c from -270^ to 
+90^ is equivalent to varying the quantity ^ - RA 0 from -180^ to +180^. 

In figure 11, the variation of tangent latitude with nodal orientation 
relative to the subsolar meridian is shown for 600-km circular orbits having 
inclinations of 50^ and 70^, with solar declination as a parameter. Fig- 
ure 11 (a) shows this variation for the 50^ orbit with solar declinations of 0 ^, 
±16^.1, and ±23^.44. At each value of - RA 0 the limiting northern and 
southern tangent latitudes are shown with 63 as a parameter. The choice of 
63 = ±16^.1 is not arbitrary; it is the limiting declination for which periods 
of continuous sunlight are possible with this orbit. For I 63 I > 16^.1, there 
are orbit plane orientations which produce no occultations. It is clear from 
figure 11 (a) that limiting declinations of ±16^.1 produce the maximum limits on 
tangent latitudes for this orbit. The fact that the maximum magnitude of the 
tangent latitude occurs at some time other than when the terminator passes 
through the poles may appear contradictory, as the equinoxes mark the only 
times at which pole-to-pole coverage is possible under any orbit conditions. 
Qualitatively, this result may be explained by noting that the poles are never 
visible from the 50^, 600-km orbit, and occultations can only occur at visible 
latitudes where the spacecraft can view the sun perpendicularly across the 
terminator. 

In figure 11(b), available tangent latitudes are shown for a 70^, 600-km 
orbit. Here, continuous sunlight can occur at any solar declination for some 
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angular orientation of the orbit plane. Also, polar coverage is available from 
this orbit when the terminator goes through the poles (i.e., when §0 = 0°). 

The solar declinations of ±10° have no special significance for this orbit and 
are just intermediate values to illustrate the tangent latitude behavior . 

As an example, consider the envelope of available tangent latitudes for 
the 50°, 600-km orbit as a function of solar declination. Minimum and maximum 
values of Xtp are available from figure 11, for the orbits and solar declina- 
tions considered there, with equation (6) available to fill in additional 
points as needed. In principle, the analytic derivative of equation (6) with 
respect to c or some trigonometric function of c could be set equal to zero 
and evaluated with the required orbital parameters, but such an analytic solu- 
tion has not been found. The maximum magnitude of tangent latitude IXipJ as 
a function of solar declination is illustrated in figure 12 for the 50°, 600-km 
orbit. One boundary on the latitude is formed by the solar terminator itself, 
while the other is formed by orbital constraints. Note that the maximum lati- 
tude magnitude of 73°. 9 is reached when j S©! = 16°.l, and this occurrence 
also marks the boundaries between which periods of continuous sunlight are pos- 
sible. As indicated on the figure, ] 6g| ^ 16°.l between May and August and 
again between November and February. 

Figure 13 shows the envelopes of available tangent latitudes for several 
600-km circular orbits which result from consideration of solar declination and 
of figure 11 in the case of the 50° and 70° orbits and similar calculations for 
the other orbits. The abscissa is linear with time, rather than with solar 
declination, so the resulting envelope could be superimposed on plots of actual 
or simulated tangent latitude variations such as figure 4(a). Table 6 (adapted 
from ref. 3) gives a tabulation of solar declination and right ascension as a 
function of time from January 1 and the vernal equinox, from which the termina- 
tor restrictions can be obtained for figure 13. Clearly, all the points north 
of the terminator curve in the Northern Hemisphere or south of terminator curve 
in the Southern Hemisphere are inaccessible with any orbit orientation. Any 
additional constraints on the envelope are due to the orbit. As previously 
noted for the 50°, 600-km orbit, the maximum northern or southern tangent lati- 
tudes occur not at the equinoxes, when the terminator goes through the poles, 
but between the times of the equinoxes and solstices. The maximum available 
latitude of 73°. 9, north or south, for the 50° orbit is equal to 90° minus the 
minimum solar declination which allows some periods of continuous sunlight. 

Note that in figures 12 and 13 the envelopes are not smooth curves; the orbit 
constraint curve intersects the solar declination constraint rather than blend- 
ing in with it. 

Also shown in figure 13 is the envelope for a 60°, 600-km orbit, which has 
a maximum available tangent latitude of 83°. 9, and a 40°, 600-km orbit which 
should never reach the solar declination limit according to figure 9. The 
envelope for a 70°, 600-km orbit follows the solar terminator curve during each 
hemisphere's summer season. This is to be expected from figure 9, which indi- 
cates that at 600 km, any inclination above 66°. 1 (90° - 23°. 9) allows periods 
of continuous sunlight regardless of solar declination. This is equivalent to 
saying that such an orbit can always be oriented so that the Sun may be viewed 
across the terminator at either its northern or southern limits (but not both 
at the same time of year) . 


14 



Finally, the maximum tangent latitudes, north or south, which are avail- 
able from various orbits are summarized in table 7. These data are just exten- 
sions of the calculations for figure 9. Also given in the right column of 
table 7 is the minimum inclination which allows polar coverage at the equinoxes 
as a function of circular orbit altitude; these data are also based on the cal- 
culations used for figure 9 and are included here for completeness. 

It can be seen now that mission-particular solar occultation data such as 
those found in figure 4 can be generalized with the help of the preceding dis- 
cussion. The- tangent latitudes vary within boundaries which are not visible in 
figure 4, and these can be specified for any required orbit, as in figure 13, 
so that at least the limiting effects of varying launch conditions can be 
assessed. These effects include a seasonal component due to the Sun's motion 
in declination and a daily component, equivalent to varying the nodal relation- 
ship of the orbit to the Sun. Information such as that contained in figure 11 
is directly applicable to determining the details of measurement opportunities 
available during a given part of the year. If periods of continuous sunlight 
are of interest, for either their deliberate inclusion or exclusion, figures 9 
and 11, or their equivalents for particular orbits of interest, provide the 
necessary mission-planning data. The geographical distribution of measurements 
can be tailored to some extent by choosing the season and time of day of launch. 
This is especially important for short-term missions, as orbital precession is 
too slow to yield in a few days the variety of measurement opportunities avail- 
able over the space of a year. 

To illustrate how tangent latitudes vary within the boundaries defined in 
figure 13, figure 14 illustrates tangent latitudes for a 50°, 600-kra orbit. In 
figure 14(a) the initial nodal orientation relative to the Sun is the same as 
for the ncsninal mission. In 14(b) the orbit plane is rotated by 90° to produce 
sunrises and sunsets at the start of the mission instead of after a period of 
continuous sunlight. The tangent latitude envelope for this orbit, from fig- 
ure 13, is superimposed on each part of figure 14; the periods during which 
continuous sunlight is possible are also shown. As expected, the resulting 
cycles of tangent latitude for each initial condition are contained within the 
envelope, and periods of no occultation are restricted to those times of the 
year when the magnitude of the Sun's declination exceeds 16°. 1. Note the 
faster cycles of tangent latitude for the 50° orbit relative to the 70° nominal 
orbit. From reference 3, the nodal precession rate is proportional to cos i: 


-2.0646 X lO''^ cos i 

» deg/day (7) 

a3-5 

• 

Table 8 gives nodal precession rates for a variety of orbits. Note that is 
zero at an inclination of 90° and positive for angles greater than 90°. Incli^ 
nations of around 100° are in the range of Sun-synchronous orbits, for which 
needs to be 0.9856473 deg/day (see ref. 3). A conflict is clearly evident 
between frequency of coverage at some latitudes, driven by the nodal precession 
rate, and total portion of the globe covered, a function of orbit inclination 
and altitude. For example, equatorial latitudes are covered on 23 occasions 
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during 1 year from a 50°, 600-km orbit, but there is no polar coverage. To 
achieve this at the same altitude, a higher inclination is required and for the 
70° nominal mission, equatorial coverage is reduced to 1 5 occasions per year. 

A final item of parametric interest relative to figure 14 is the distribu- 
tion of measurements as a function of latitude. It is clear from figures 4(a), 
(c) , and (d) that there are many more measurements at some latitudes than at 
others. Distributions of number of measurements as a function of tangent lati- 
tude, in 5° bands, are shovm in figure 15 for 1-year missions with 600-km 
orbits having inclinations of 70°, 60°, and 50°. The peaks in the northern and 
southern latitudes, away from the equator, are typical of many kinds of Earth- 
pointing measurements made from space. The lack of complete symmetry about the 
equator is due to having picked a particular initial orbital relationship to 
the Sun, the effects of which are not overridden in a 1-year mission. The ini- 
tial nodal orientation for all three orbits in figure 15 is the same as that 
used for the nominal mission of figure 4. The data for the 70° orbit may be 
obtained directly frcan figure 7 by summing all the "boxes" in each 5° latitude 
band. The total number of measurements resulting from a 1-year mission 
increases as the inclination decreases, as the total length of all periods of 
continuous sunlight decreases with decreasing i (due to the more rapid orbit 
plane precession rate) . Of course, the total number of occultation opportuni- 
ties is limited to twice the number of orbits, a value achieved when no periods 
of continuous sunlight are possible. 

In addition to considerations of coverage on the Earth's surface, through 
the tangent points and their projections through the atmosphere, there is also 
some potential interest in what is viewed on the surface of the Sun. For 
example, what are the heliographic (Sun-centered) coordinates of the unit vec- 
tor ^ which has been used in the previous mission analyses for defining 
Sun viewing conditions? If the Sun's disk is scanned from edge to edge, what 
are the heliographic coordinates available to the viewer? This may be of con- 
cern in avoiding or accounting for sunspot activity, which is known to be a 
function of heliographic latitude (see ref. 5, for example). For instruments 
whose field of view is small compared with the angular size of the Sun (approx- 
imately 0.5 deg) it is also possible that accurate instrument calibrations and 
data interpretation may be compromised if the viewing conditions on the Sun are 
not well known during the time of calibration and data taking. Such informa- 
tion can be generated during preliminary mission analyses by making use of 
information calculated at the time of occultation. The problem is complicated 
geometrically by the Sun's rotational axis not being perpendicular to the 
ecliptic; that is, its equator does not lie in the ecliptic. Earth-centered 
vectors must first be transformed into a heliocentric-ecliptic system and then 
into the heliographic system, or vice versa. Because of the peripheral nature 
of this problem, its self-contained solution is dealt with in appendix B. An 
example is given there of how the heliographic coordinates might change over 
the year during the nominal mission. 

STATISTICAL INTERPRETATION OF DATA FROM SOLAR OCCULTATION MISSIONS 

In order to understand how data from long-term solar occultation missions 
may be interpreted, it is necessary to augment the information already generated 
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for the simulated missions described previously. The spatial and temporal 
distributions of measurement opportunities have already been specified for 
the nominal mission of figure 4 (70° inclination, 600-km altitude, 1 year) . 

In keeping with the general approach outlined in the introduction to the paper, 
a single additional "output" quantity q is now associated with each measure- 
ment opportunity. For the nominal mission there are 9480 such opportunities. 
The parameter q is a function of tangent longitude, tangent latitude, and 
time. It is assumed for the present purposes to have been obtained directly 
without measurement error, even though the process of obtaining values of q 
as a result of an actual mission may well be an arduous and error-prone task. 
The physical variable which q represents is assumed to have both systematic 
and random variability. No other physical characteristics need be attributed 
to q. A dimensionless numerical value is assigned to each sunrise or sunset 
measurement opportunity, uniquely specified by its time and position coordi- 
nates. The mathematical model which produced the simulated measurements is not 
given except as it can be determined from the measurements themselves. This 
approach is an attempt to free the discussion which follows from the prejudice 
of additional knowledge which may not exist for a real mission. 

The problem of data interpretation is to group the measurements in a mean- 
ingful way and to search these groupings for patterns which emerge. It is of 
interest here to produce a function which shows how q, or its average, varies 
as a function of time. Figure 4(a) shows that the tangent point passes through 
a particular latitude several times per year, depending on the latitude. This 
occurs on 15 different occasions per year at equatorial latitudes (for the nom- 
inal mission), but only once or twice at polar latitudes. Additionally, fig- 
ure 6 has shown that complete longitudinal coverage is obtained at about 25° 
intervals, during a period of time (a little over a day) in which the tangent 
latitude changes very little. Thus, figures 4(a) and 6 can be used to justify 
grouping solar occultation data in subsets corresponding to measurements taken 
during small time intervals, within bands of latitude. Longitudinal dependence 
may be removed if it can be assumed to be negligible, or to be properly aver- 
aged by considering the available data within a band of latitude. In fact, for 
many measurements of interest, dispersion and transport of atmospheric constit- 
uents is much more rapid to the east and west than it is to the north and 
south, leading to longitudinal uniformity when no comparable uniformity 
exists as a function of_latitude. In any event, it is necessary to quantify 
this assumption. Let t be an average time during which a subset of mea- 
surements is made, and let X be the midpoint in a latitude band (5° bands 
are used for this study) ; 


t = 


^last ^first 
2 


^high ^low 
2 


( 8 ) 
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within each latitude band and time interval, there exists a true mean Pq; 


Mq(X,t) 


1 1 1 p p 

. . = — T- 77 \ q(X,L,t) dX 

X,t 27T At AX Jq ^At ‘^X 


dt dL 


( 9 ) 


It is assumed that for ng measurements of some parameter q within a partic- 
ular latitude band, the true mean can be approximated as a function of latitude 
and time by 


q(X,t) 


X,t 


nq 

5" qi(X,L,t) 
i=l 


nq 


( 10 ) 


Table 9 gives a set of values of q(X,t) for the nominal mission. Within 
each 5° latitude band, starting at the South Pole, measurements are grouped 
according to time in the following way: if a measurement time differs from 

the preceding one by less than a day, it is retained in the current data sub- 
set; otherwise a new subset is started. This procedure groups the data in the 
natural way indicated by figure 4(a). Near the equator there are 15 such sub- 
sets in a year. Table 9 lists for each subset (there are from 1 to 1 5 subsets 
in each latitude band) the number of measurements nq, the mean value q, the 
sample standard deviation Sg, the average time of the measurement t, and the 
first and last measurement time ^first ^last referenced to the begin- 

ning of a calendar year, January 1, 1981. It has been assumed for simplicity 
of presentation that the data are periodic over a year so that information 
obtained on any date is the same except for random components as that obtained 
on the same date in adjacent years. Hence the data obtained in 1982 can be 
used to fill in the 1981 data prior to the mission launch time. This need not 
apply to an actual mission, but no difficulty is introduced into the analysis 
if such an assumption is not made. The sample variance Sq^ is computed in 
the usual way: 



^ (qi - 5)^ 

i=l 


nq - 


1 



( 11 ) 


For those cases where nq = 1 , Sq is not defined, and such places are 
noted by blanks in table 9; the first and last measurement times are also left 
blank. 
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There are still uncertainties to be resolved about each q and its rela- 
tionship to the true mean Pq over the same time interval. In general, the 
distribution of various atmospheric constituents represented here by q may be 
driven by diurnal and short-term cycles in addition to a longer term (seasonal) 
cycle whose behavior is sought as a result of the proposed 1-year nominal mis- 
sion. It is clear that diurnal cycles are fundamentally inaccessible through 
solar occultation measurements, which are restricted to local dawn or dusk con- 
ditions. The short-term cycles might be driven in longitude by what are com- 
monly known as planetary waves (which may also exhibit long-term patterns) , and 
their behavior might typically be of considerable interest in an actual mission. 
The sampling frequency of q (about 30 per day) or, equivalently, the longitu- 
dinal spatial resolution of about 25° (recall fig. 6) poses certain problems 
for extracting short-term temporal information from this type of mission. 
Detailed consideration of these matters is beyond the scope of this paper. 

Note that the span of time over which each q is formed, t]^ast “ ^first' 
short compared with the seasonal cycle of interest and it is assumed that the 
behavior of the short-term cycles is such that q is equal to the true mean 
as defined in equation (9) . If this assumption breaks down because of planetary 
wave behavior, other short-term behavior, or more permanent longitudinal nonuni- 
formities, then a one-dimensional spatial concept of latitude zonal averages 
over the time scales of interest would not be valid and should be discarded in 
favor of a two-dimensional spatial grid. These questions deserve closer scru- 
tiny when measurements are considered for atmospheric constituents whose actual 
behavior is known to some extent, either from theory or previous measurements. 

In particular, the restriction of solar occultation measurements to dawn and 
dusk means that some understanding of diurnal behavior is essential for proper 
interpretation of the q's. 

Suppose now that it is desired to determine a yearly mean Pq as a func- 
tion of latitude: 


Pq(>^) 


1 

1 1 

p2n 

p365 r 

X 2tt 

365 AX * 

1 ' 

\ \ 


'AX 


( 12 ) 


It is certain that pg cannot be obtained just by taking the unweighted mean 
of all the q's in each latitude band. The difficulty with such a procedure 
is easily seen by referring to the column labeled nq in table 9. For example, 
in the latitude band from -65° .to -60°, nq varies from 5 to 218, so that an 
unweighted average of all the q's in this band favors some measurement oppor- 
tunities over others when there is no reason to do so. This difficulty is 
found in all the latitude bands, so taking unweighted averages is clearly not 
an acceptable use of the available data. Another approach is to use the q's 
themselves as data and average them to estimate pq. This could be thought 
of as a way of weighting the individual points within each subset of measure- 
ment opportunities. It is an intuitively attractive way to achieve the desired 
results, especially when there appear to be "enough" values of q. If the 
q's are used in this way, it is then tempting to say, as in equation (10), 
that for each latitude band the yearly average Pq is estimated by a 
quantity Q: 
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(13) 


Q(X) 


X 


nq 


nq 


where ng is the number of data subsets for each latitude band in table 9. If 
this is an acceptable procedure, then the variance is 


nq 

ZIqi2- 

i=l 


Hq — 1 



However, it is a fundamental characteristic of the solar occultation tech- 
nique that such calculations are not justified in the statistical sense. 
Equations (13) and (14) presume that the samples (the q*s) are taken uni- 
formly (in time) from the total population. The distribution of q appears 
uniform near the equator but is clearly not uniform near the poles. At inter- 
mediate latitudes the distribution is difficult to assess in a strictly quanti- 
tative way. Because of the problems with applying otherwise straightforward 
descr iptive^statistics, no statistically acceptable confidence limits can be 
placed on Q as calculated with equation (13). That is, it is not necessar- 
ily true that there is a 95-percent probability that the true mean yg differs 
from Q by no more than ±2sg, as would be the case if the q*s were normally 
distributed. The well-known t-statistic, which might be used for establishing 
confidence limits when the population variance is not known, is not applicable 
either, as it, too, requires that the samples from the population be randomly 
selected. 

An alternative way to arrive at the desired goal of establishing a yearly 
mean within a latitude band is to propose a model for describing the behavior 
of q. At the center of each latitude band let 


q = + A sin (360t + y) + e (15) 

where ^ is the mean value of q averaged over a year, A is an amplitude, 
t is time in fractions of a year, y is a phase factor (in degrees), and e 
is a residual term. This is a very simple periodic model which, through y, 
allows some adjustment of the location of the maximum and minimum values within 
the yearly cycle. The residual term can incorporate random error as well as 
unmodeled variability. Equation (15) can be rewritten as follows: 
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q = Q. + A cos Y sin 360t + A sin y cos 360t + G 


(16) 


Or, letting Ci = C2 = A cos y, and C3 = A sin y, 


q = C-| + C2 sin 360t + C3 cos 360t + e 


(17) 


For n measurements of q, M-| , . . 


Ml 




M, 


n 


1 sin 360ti cos 360ti 


1 sin 360ti cos 360t£ 


1 sin 360tjj cos 360tjj 


C] 

C2 

\1^3J 


(18) 


In matrix notation. 


M = ^ (19) 

where C is a 3 x i matrix of parameters to be estimated and B is a 3 x n 
matrix of coefficients evaluated at the measurement times t-| , . . .,ti,. . .,tn. 
It is then well known that the set of parameters which minimizes the sum of the 
squares of the residuals (the unweighted least squares estimate) is given by 


8 = (bTb)“''bTm 


( 20 ) 


✓s 

From the elements of C (C-| , 


C2 , and 


C3) the desired model parameters are 


^ = Cl 


A 


Y = tan""* (C3/C2)> 
A = C2/C0S y 


(21) 
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Recall that the problem with calculating statistical means, as in equa- 
tion (13), was one of interpreting the results. Now, with a model and its 
parameters, this interpretation is made straightforward in its own restricted 
way. In the least squares sense described above, the model parameters in 
table 10 are the best ones obtainable for the model used. It remains only to 
calculate how much of the observed variability in the data (the q's) is 
represented by the model. With three parameters, three data points can be 
represented perfectly, while more than three are generally imperfectly repre- 
sented. This may be true even if the postulated model is the true model, 
because of the possible existence of random error in the data. It can be shown 
that for any observed value from a set of observed values, its expected 

value Yi (as calculated from a model) , and the mean value Y of all 
observed values in the set, the following relation holds: 


(Yi - Y)2 = XI (^i - Y)2 + XI (^i - ^i)^ (22) 

i i ^ . 

Sum of squares = Sum of squares due + Sum of squares about 
about the mean to regression the regression 


This well-known statistical result expresses the fact that the variability 
about the mean (the sum of squares about the mean) is explainable in terms 
of the modeled variability (the sum of squares due to regression) plus the 
unmodeled variability (the sum of squares about the regression) . It is 
desired that the unmodeled variability be small, that is, that the ratio 



1 


R (23) 

XI (Yi - y)2 


be close to unity. If R is 1 (or 100 percent), then the model perfectly 
represents the data. If R is less than 1, then there is variability which 
has not been accounted for with the existing model; this may be random or 
systematic variability. 

The last column of table 10 shows R for the least squares fit to each 
applicable latitude band. It is interesting to note that almost none of the 
variability in the data is represented by the model near the equator. This 
corresponds to the decreasing amplitude of the sine curve. A, near the 
equator; there is very little systematic variability of the form postulated. 

It is tempting to conclude that the variability which is seen near the equator 
is only random variation in q, but this conclusion, based on the information 
presented here, is unwarranted; it can be said only that the postulated model 
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does not work everywhere. Farther away from the equator, the model appears to 
do better, indicating the possible existence at these latitudes of a marked 
annual cycle of simple (sinusoidal) structure. This model seems at least con- 
ceptually to have obtained a result, the yearly average value which was not 
obtainable with the simple descriptive statistics described previously. It 
could be argued that even better results could be obtained with a more sophisti- 
cated model. However, there are some basic problems with the present approach 
which are analogous to the lack of confidence in computed statistical means. 

The least squares technique guarantees, in a particular sense, the best values 
for the model coefficients, given a particular set of data. But, the goal of 
modeling is not just to represent one set of data, but to represent the func- 
tional behavior of the quantity of interest. If the same model were applied to 
other samples of q, different coefficients would result. In fact, if sev- 
eral sets of q's are used, all of which are uniformly distributed over a 
year within each latitude band, the resulting coefficients can be shown to be 
normally distributed with variances related (through a measure of the error 
variance) to the diagonal components of the matrix (B^B)~^ as defined in 
equation (20) . These variances are largest when ng is smallest. When 
ng = 3, the data are represented perfectly by the model, but the variances 
associated with the coefficients are no longer defined. This means, for exam- 
ple, that the behavior of q cannot be predicted anywhere else within that 
latitude band. Closer to the equator, the variances are smaller, and hence the 
model has more usefulness as an interpolative tool. The coefficients can be 
assumed to be correct and the residual term can be assumed to represent random 
components with variance where s^^ is the sum of squares of the resid- 

uals based on three degrees of freedom (because three parameters have been 
calculated) : 


21 (qi - 5i)2 

i=1 

Sr2 = ^ (24) 

nq - 3 


The estimated value of is obtained from equation (15) using the coef- 

ficients from equation (21). Despite the difficulty of using the model coeffi- 
cients with Sj^^ to predict a range for the values of q, sj^^ is given in 
table 10 for each latitude band having enough data; the values tend to be larg- 
est near the equator, but there is no explanation for this behavior within the 
information presented. 

There are two reasons for not calculating the coefficient variances in 
this paper. First, the samples of q are not demonstrably uniformly dis- 
tributed in all cases, especially near the poles. Second, the error variance 
associated with the q's is not known. It may be related to Sj^^^ which 
would be the variance due to random components in q if the postulated model 
were the true model and the residual represented only random contributions to 
q due to the properties of the physical parameter represented by q. (It 
may be believed that the distribution of any quantity in the atmosphere is 
ultimately deterministic, but the time scale of the driving processes may be 
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such that the results would appear truly random when superimposed on the 
occultation measurements.) It is more likely that the model is too simple, 
incomplete, or wrong, and that the error variance should be some composite of 
variances of the q's (the original data) as determined from equation (11). 
That is, there should be some global function which determines the size of the 
random components superimposed on q. A measure of the true error variance 
is rarely possible; it requires repetition of measurements in both space and 
time. This is fundamentally unattainable with satellite measurements and 
requires, as a minimum, the assumption that equivalent future times can be 
found to substitute for repetitive measurements at the same time. 

As a final comment on the use of models to represent satellite data, note 
that the success of such models in representing data no more guarantees their 
physical correctness than a lack of such success indicates their unsuitability. 
In the present example, data have been represented very well over part of the 
globe (near the poles) , but the data themselves are insufficient for qualita- 
tively satisfying explanations of parametric behavior. In other regions, 
common sense seems to dictate that explanations in terms of small systematic 
variability and large random components are reasonable in spite of a lack of 
statistical rigor in the arguments. Further definition of this situation 
demands additional insight from theory and experiment as well as ample justifi- 
cation for detailed analysis. 


CONCLUDING REMARKS 

In the first part of this paper, the orbit dynamics of solar occultation 
have been discussed first for a particular mission and then in a more general 
way. This order of presentation uses a hypothetical mission to define impor- 
tant characteristics of all solar occultation missions so that the most signif- 
icant points can be amplified in the parametric discussions. Observability of 
the Sun from space as it rises or sets on the horizon is determined by orbit 
parameters and seasonal factors; the interdependence between satellite and Sun 
is a complicated one, as both source and observer contribute to a constantly 
changing geometry. Emphasis has been placed on understanding the distribution 
of measurement opportunities over the surface of the Earth. The critical spa- 
tial coordinate is latitude. Good longitude coverage is guaranteed by the 
Earth's rotation, at a much higher rate (1 revolution per day) than the rota- 
tion of the orbit plane in inertial space (a few degrees per day) . Pole-to- 
pole latitude coverage is available only in the spring and fall, when the solar 
terminator passes through the poles. Even then, orbit parameters - altitude, 
inclination, and orbit plane orientation - may prohibit coverage at the north- 
ern and southern latitude extremes. A year-long nominal mission (70° inclina- 
tion, 600-km altitude) has been chosen to provide coverage which extends into 
both polar regions in the spring and fall. Solar occultation latitude coverage 
patterns are characterized by their cyclic behavior, driven by precession of 
the orbit plane and the Sun's motion in inertial space. At some higher incli- 
nations and latitudes, periods of no occultation can occur, during which the 
Sun is always in sight of the spacecraft. These periods typically last a few 
days. A special case occurs for Sun-synchronous orbits, where the precession 
rate matches the average angular rate of the Sun (about 1 deg/day) . Then the 
initial geometry is maintained, on the average, and geographical coverage is 


24 



limited in latitude or nonexistent if the geometry is such that no sunrises or 
sunsets may be viewed. 


Temporal coverage is characterized most prominently by its restriction to 
local dawn or dusk. With 2 measurement opportunities per orbit, there can be 
on the order of 10^ measurements during a year. The uneven seasonal distribu- 
tion of the measurements (within bands of latitude, for example) is such that 
straightforward statistical interpretation of time-averaged quantities of 
interest is not possible. This problem has led to the second major division 
of the paper, in which some statistical limitations of solar occultation data 
are discussed in detail. The procedure involves the association of a single 
output quantity with each measurement opportunity. All that is known about 
this quantity, for the purposes of this paper, is that substantial variability 
is revealed by the available data. Although there are certain qualitatively 
"obvious" conclusions which can be drawn from an examination of the data, the 
nature of the data prevents statistically meaningful confidence limits frcxn 
being placed on some apparently simple calculations. If descriptive statistics 
are discarded in favor of mathematically modeling the functional behavior of 
the observed quantity, it is shown that the extent to which variability of a 
particular finite set of data can be represented in a mathematical way has 
little to do with the physical appropriateness of the model. (The simple model 
proposed represents the data quite well over some portions of the globe, but 
there is no reason whatsoever to embrace its physical implications.) If it 
is realized that functional representation of certain behavior is the goal of 
modeling, and not just the representation of a specific set of data, it has 
been concluded that statistical confidence in the modeling coefficients may be 
as lacking as confidence limits were for averaged quantities in the descriptive 
statistics. The persistent interpretive problems with these data are reassur- 
ing in a way, because they arise from insufficient information about the quan- 
tity of interest; any technique would be suspect if it appeared to overcome 
such difficulties in the absence of additional information, either actual or 
assumed. 

This paper shows that long-term solar occultation missions can contribute 
significantly to understanding the Earth's atmosphere. In some cases, when 
measurements obtained in this way will be the only ones available, their sta- 
tistical adequacy is of secondary importance. For valid statistical analyses 
to be applied to global descriptions of parametric behavior, such measurements 
must eventually be supplemented with other sources of information. Similar 
approaches can be applied to occultation involving objects other than the Sun, 
natural or man-made. The caveats encountered in data analysis are applicable, 
as well, to other environmental measurements; they point to the need for care- 
ful experiment design whenever such control can be exercised. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
January 10, 1980 
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APPENDIX A 


DETERMINING THE CONDITIONS UNDER WHICH THE SUN IS OCCULTED BY THE EARTH 

RELATIVE TO A SATELLITE HAVING FIXED ORBIT ELEMENTS 

Solution of the occultation problem involves determining the points on 
its orbit at which a satellite enters or leaves the Earth's shadow. The 
geometry is shown in sketch (A1 ) . The symbols used for this appendix form a 



y 


Sketch (Al) 

self-contained set, and they are not all defined in or consistent with the 
symbol lists in the body of the paper or in appendix B. All the vector quan- 
tities are expressed in a system whose x-y plane lies in the plane of the sat- 
ellite orbit, with x and y oriented as indicated in sketch (Al ) . Such 
a coordinate system is called the PQW system. Vectors can be rotated into the 
PQW system by standard procedures. (See, for example, ref. 6 .) 

Aq to A4 constants as defined below 
->■ 

d vector perpendicular to the terminator, extending from the terminator 

to the satellite, km 

e eccentricity 
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f true anomaly, deg 

p semilatus rectum, km 

->• 

r radius vector to satellite, km 

?0 unit vector to the Sun 

rg Earth radius vector to the terminator, km 

X unit vector in direction of the perigee of the satellite orbit 

y unit vector at right angles to x along the semilatus rectum 

il^ angle between fg and r, deg 

The geometric constraint is the requirement that at a point where the satellite 
enters or leaves the Earth's shadow, 

fg • d = -d = “(r^ - r0^)V2 = j- cos ^ 

That is, T0 and d are antiparallel. This restriction can be written in 
another way; 

?0 • d = r© ♦ (r - r0) = ?0 • r = r?0 • {x cos f + y sin f) = r cos 4) 


The X- and y-components of the unit vector to the Sun are 


Then, 


/N /N 


r©x = r© ♦ X 


/N 


r©y = r© • y 


r©x cos f + r©y sin f = - 


(r2 - r ® 2)'/2 


Since r = p/(l + e cos f) , 


1 + e cos f 


(r©x cos f + r©y sin f) + 


lV2 


(1 + e cos f ) 2 


- r®" 


= 0 
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This can be rewritten in terms of cos f; 


Aq cos^ f + A-| cos^ f + A2 cos^ f + A3 cos f + A4 = 0 


The constants are 



A2 = 






2 2 2 2 2 
+ 2(r0y - r0x) (1 - t0y) - 4r0xt0y 




A4 



2 2 
>^0y) + (1 ~ >^0y) 


2 


The fourth-order equation may be solved by standard procedures. The two 
spurious roots may be rejected by noting that |r0 x ?’| = r sin = I?©! and 
r© • r/r ^ 0 for occultation to occur. 
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VIEWING POINTS ON THE SUN'S SURFACE FROM A SPACECRAFT 

The purpose of this appendix is to develop the equations necessary for 
locating a point on the Sun's surface as it is viewed from a spacecraft in 
orbit around the Earth. The symbols used for this appendix form a self- 
contained set, and they are not all defined in or consistent with the symbol 
lists in the body of the paper or in appendix A. 

a central angle between the intersection of the scanning plane with the 

solar equator and the sub-Earth meridian, deg (see sketch (B7) ) 

A constant, defined in equation (B7) 

B heliographic latitude of the viewing point, deg 

Bq heliographic latitude of the sub-Earth point, deg 

c central angle between the intersection of the scanning plane with the 

solar equator and the sub-Earth point, deg (see sketch (B7) ) 

ip angle between the scanning plane and the solar equatorial plane, deg 

ig inclination of the solar equator to the ecliptic, deg 

J.D. Julian date 

difference in heliographic longitude between the meridian passing 
through the viewing point and the intersection of the scanning 
plane with the solar equator, deg (see sketch (B7)) 

L heliographic longitude of the viewing point, deg 

Lq heliographic longitude of the sub-Earth point, deg 

M heliographic longitude of the ascending node of the solar equator on 

the ecliptic, deg 

✓X 

P unit vector in the direction of the Sun's north pole, may be specified 

in an ecliptic coordinate system (P^c) 3 right-ascension— 

declination (Earth equatorial) system (Peq) 

->■ ^ 

r(r) position vector from Earth's center (or unit vector) of a 

spacecraft, km 

rrp vector from Earth's center to a point of tangency at the Earth's 

surface or in its atmosphere, km (see sketch (Bl)) 


✓X 

r fp ( r / r ^ r ) 


unit vector in the direction of rci 
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S sub-Earth point on the Sun 

S vector from the spacecraft to the point of tangency 

SC spacecraft 

V viewing point on the Sun 

X0(X0,y0,Z0) unit vector to the Sun 
X matrix defined in equation (B3) 

a angle between S and -r (see sketch (B1 ) ) 

e obliquity of the Earth's equator to the ecliptic, deg 

n(Tix/ny,nz) unit vector, rq. x x© 

rip unit vector, p x x© 

0 angular part of a polar coordinate system for describing viewing 

points on the Sun's surface (see sketch (B2) ) 

A heliocentric longitude of the Earth, deg 

p number between 0 and 1 which locates the viewing point relative 

to the center of the solar disk (see sketch (B6) ) 

p(Px,Py,P2) unit vector defining the orientation of the scanning plane 

relative to the spacecraft 

T direction of vernal equinox vector 

central angle, sin"l P 

celestial longitude of the ascending node of the solar equator on 
the ecliptic, deg 


Subscripts; 

x,y,z components of a unit vector along Earth equatorial right-ascension— 
declination coordinate system 

Assume that the Sun is viewed frcan a spacecraft along a line of sight 
directed at the center of the solar disk. The "cross hairs" of this viewing 
system are the projections of r.ji and fj on the solar disk, where the 
geometry is as follows: 
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Sketch (B1 ) 

Points on the solar disk are located by the polar coordinate system (p,9) as 
shown below: 



The unit vector p in the right- ascension— declination system is found by 
solving the following set of equations, where 9 is specified: 

/N _ 

p • X0 = 0 

p • = cos 9 > (Bl) 

p * n = -sin 9 
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Or, in matrix form: 


where 


px = 

[0, 

cos 9 , 

-sin 0] 

(B2) 



*^Tx 

Px 



X = 

y© 

rTy 

ny 


(B3) 


_z© 

^Tz 

Pz^ 






Then, 

p = [0,cos 9, -sin 9]x“"* 
Algebraically, , 

Px^0 + PyyQ + PzZQ = 0 
Px^Tx PyJ^Ty Pz^Tz ~ cos 0 ^ 
Pjjflx + Py0y + Pz’^Z = ^ 

and the components of p are 


Pjj = — 1^- (r^ sin 9 + Py cos 9)z0 - (rip 2 sin 0 + ri^ cos 0)y0j^ 


Py = ~[+(tTx sin 9 + Tijj cos 0)Z0 - (rijig sin 9 + ^2 cos 9)X0j ^ 
P 2 = ~j^“(tTx ® ■'■ 9x c°® ®)y© “ (J^Ty 9 + Py cos 9)X0j 


(B4) 


(B5) 


(B6) 
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where 


A = (n^r^y - nyrTx)Z0 + (’^z'^Tx * Vtx)^© + (Vtz " VTy)^ (^7) 

The Sun’s coordinate system is detailed in reference 7. Specifically/ the 
inclination of the solar equator to the ecliptic and the celestial longitude of 
the ascending node of the solar equator on the ecliptic are 


i0 = 7°. 25 (B8) 

^0 = 73°. 666667 + 3.821583 x 10"5(J.D. - 239 6759.5) (B9) 


where 239 6759.5 is the Julian date (J.D.) for 1850.0 (see ref. 3 for an 
explanation of time and coordinate systems) . A heliographic longitude- 
latitude system is established on the Sun's surface as follows: 



Sketch (B3) 

The heliographic longitude of the solar equator ascending node on the ecliptic 
M is arbitrarily referenced from a solar meridian that passed through the 
solar equator ascending node at Greenwich mean noon, January 1, 1854, 

J.D. 239 8220.0. The sidereal rotation period of 25.38 mean solar days is 
used to establish the Sun's rotation rate of 14.18439716 deg/day, and M 
is given by 

M = 292°. 766 + 14.18439716(243 0000.5 - J.D.) (BIO) 
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The unit vector for the orientation of the Sun's pole in ecliptic coordinates 
^ec 


^eCfX “ sin i 0 sin 

Pgcy = “Sin i© cos % / (Bll ) 

^ec,z “ ^0 > 


The central meridian is that meridian which passes through the sub-Earth point 
S. For all practical purposes S is identical to the subsatellite point for 
any near-Earth orbit. The sub-Earth coordinates are calculated from the 
Earth's heliocentric longitude A: 



sin(S ^0 - A) sin i© = sin Bq (B 12 ) 

cos{Qq - A) = cos (90° - Bq) cos 90° + sin (90° - Bq) sin 90° cos (M - Lq) 

cos(fl© - A) = cos Bq cos (M - Lq) (B13) 

sin (90° - Bq) sin (f2@ - A) 


sin (90° - i© sin (M - Lq) 


cos Bq sin (M - Lq) = sin (fl© - A) cos i© 


(B14) 


34 



APPENDIX B 


Equations (B12), (B13), and (B14) are sufficient to solve unambiguously for the 
coordinates Lq and Bq. 

The projection of P across the solar disk defines a plane which cuts 
the Sun in a great circle through the sub-Earth point. The geometry is shown 
below: 



The viewing instrument may scan along the projection of p to some point V. 
The angle which the scanning plane makes with the solar equator is required. 

In an Earth equatorial right-ascension— declination system, this plane is char- 
acterized by the unit vector fip: 


Tip = P X X0 


(B15) 


Then, the angle between the equator and the scanning plane ip is 


ip = cos ■> (fip • Pgq) 


(B16) 


where the unit vector for the Sun's north pole has been transformed from an 
ecliptic system (see eq. (Bll)) into the Earth equatorial right-ascension— 
declination system: 


Peq,x = sin i© sin ^ 

Peq,y = -sin i© cos cos e - cos i© sin e > (B17) 

Peq,z “ -sin i© cos sin e + cos i© cos e 

J 


The Earth's obliquity £ is about 23°. 44. 
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To define the viewing point V along consider moving a certain dis- 
tance away from the center of the solar disk. The geometry is shown below; 



Sketch (B6) 

Some fraction of the solar disk, given by p between 0 and 1 , is converted to 
a central angle on the surface of the Sun through the relationship for 


The problem is solved now, in principle. Typically, the solution is coded 
for computer implementation and it is necessary to exercise caution in applying 
spherical trigonometric relationships. The geometrical relationships are sum- 
marized below: 
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Interpretation of these angular quantities is influenced by three factors; the 
angle between rip and P; the angle between r and P (where r, rip, and ^ 
must be expressed in identical coordinate systems, the Earth equatorial right- 
ascension— declination system, for example) ; and the value of sub-Earth latitude 
(specifically, whether Bq is positive or negative) . Table B1 summarizes the 
possibilities. The following trigonometric and logical relationships may be 
applied: 


ip = cos"! (rip • &) 


i J‘5 

^ ^180® - ip 

(ip ^ 90®) 
(ip > 90®) 

c = sin“^ (sin Bo/sin ip) 


a = sin“^ (tan Bg/tan ip) 


B* = sin”^ |sin ip sin + sign ( 

r • ^)c]| 

L* = sin“^ (tan B*/tan ip) 


B = sign (r • ^)B* 



(B18) 


L = Lq + sign (cos ip)L* - sign (r • §)a 




where the sign ( ) function is equal to +1 or -1 depending on the sign of the 
argument inside the parentheses. 

An application of the set of equations (B18) to the nominal 70®, 600-km 
solar occultation mission described in the main body of the paper (see fig. 4) 
is illustrated in figure B1 . Here, the heliographic latitude for viewing the 
"top" edge of the Sun's image (6=0, P = 1) and the sub-Earth point are 
plotted as a function of time from launch for 1 year. As an example of their 
use, such data, or tabulated versions of them, could perhaps be the basis for 
a statistical analysis of the effects of sunspot activity on measurement 
opportunities. 
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0 30 60 90 120 150 180 210 240 270 300 330 360 

Time from launch, days 


Figure B1 Heliographic latitude of the subsatellite point and the top edge 
of the solar disk as viewed during sunrises and sunsets on a 1-year solar 
occultation mission using a 70°, 600-km circular orbit. 


39 



REFERENCES 


1. Harrison, Edwin F.; Green, Richard N.; Brooks, David R.; Lawrence, 

George F.; and McCormick, M. Patrick: Mission Analysis for 

Satellite Measurements of Stratospheric Constituents by Solar 
Occultation. Paper presented at the AIAA 13th Aerospace Sciences 
Meeting (Pasadena, California), Jan. 1975. (Available as AIAA 
Paper 75-57.) 

2. Russell, James M., Ill; Park, Jae H.; and Drayson, S. Roland: Global 

Monitoring of Stratospheric Halogen Compounds From a Satellite Using 
Gas Filter Spectroscopy in the Solar Occultation Mode. Appl. Opt., 
vol. 16, no. 3, Mar. 1977, pp. 607-612. 

3. Brooks, David R.: An Introduction to Orbit Dynamics and Its Application 

to Satellite-Based Earth Monitoring Missions. NASA RP-1009, 1977. 

4. Melbourne, William G.; Mulholland, J. Derral; Sjogren, William L.; and 

Sturms, Francis M. , Jr.: Constants and Related Information for Astro- 

dynamic Calculations, 1968. Tech. Rep. No. 32-1306 (Contract NAS 7-100), 
Jet Propulsion Lab., California Inst. Technol. , July 15, 1968. 

5. Allen, C. W. ; Astrophysical Quantities. Second ed. , Athlone Press, c.1963. 

(Available from Oxford Univ. Press, Inc., New York.) 

6. Escobal, Pedro Ramon: Methods of Orbit Determination. John Wiley & Sons, 

Inc., c.1965. 

7. Explanatory Supplement to the Astronomical Ephemeris and the American 

Ephemeris and Nautical Almanac. H.M. Naut. Aim. Off., 1961. 


40 



TABLE 1 ORBIT PARAMETERS FOR NOMINAL SOLAR OCCULTATION MISSION 


Semlmajor axis, a, km 6978 

Eccentricity, e 0 

Inclination, i, deg 70 

Period,^ sec ^ 5805 

Nodal precession rate, fl, deg/day -2.4868 


%odal period from first-order perturbation theory; 5805-sec period 
results in 14.883 orbits/day. 


TABLE 2.- SLOW PASSAGES OF SUN THROUGH ATMOSPHERE NEAR BEGINNING OF 


NOMINAL SOLAR OCCULTATION MISSION 


Orbit 

Minimum solar 

Time below 

number 

altitude/ km 

100 km,^ min 

102 

52.4 

9.1 

103 

42.0 

9.8 

104 

31 .6 

10.5 

105 

21 .1 

11 .2 

106 

10.5 

11.9 

107 

-.2 

12.6 

108 

-11.0 

13.3 


®The center of the solar disk is at an altitude of less than 100 km 
above the local horizon. 
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TABLE 3.- TYPICAL LIMITING CASES FOR VIEWING THE SUN IN THE VICINITY OF 
THE HORIZON, FROM THE NOMINAL MISSION 

Slow® Fast^ 

Time for measurement while center of Sun is less than 

100 km above local horizon, min 12.5 0.81 

Change in surface coordinates of tangent point: 

Longitude, deg 18 0.17 

Latitude, deg 45 <0.01 

Change in pointing angles to center of Sun: 

a, deg 2 2 

3, deg ............... 20 <0.01 

^Center of Sun moves from altitude of 100 km to the horizon and back up 
to 100 km. 

^Sunrise or sunset when Sun's vertical motion near the horizon is 
maximum. 
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TABLE 4.- DISTRIBUTION OF APPARENT VELOCITY OF SUN'S IMAGE RELATIVE TO HORIZON 


AS VIEWED FROM A SPACECRAFT DURING NOMINAL SOLAR OCCULTATION MISSION 


1 Vrell r 
km/ sec 

Number of 
events 

Cumulative 
percentage 
of total 
events 

Measurement 

time,^ 

sec 

0-0.1 

6 

0.06 

00-1 000 

0.1-0. 3 

102 

1 .1 

1000-333 

0.3-0. 5 

152 

2.7 

333-200 

0.5-0. 7 

170 

4.5 

200-143 

0.7-0. 9 

202 

6.7 

143-111 

0.9-1 .1 

224 

9.0 

111-91 

1 .1-1.3 

248 

11.6 

91-77 

1 .3-1 .5 

502 

16.9 

77-67 

1 .5-1 .7 

430 

21 .5 

67-59 

1 .7-1 .9 

844 

30.4 

59-53 

1 .9-2.1 

742 

38.2 

53-48 

2. 1-2. 3 

714 

45.7 

48-43 

2. 3-2. 5 

770 

53.9 

43-40 

2. 5-2. 7 

892 

63.3 

40-37 

2. 7-2. 9 

1156 

75.5 

37-34 

2. 9-3.1 

2326 

. .. 1 

100.0 

34-32 

Total . . . 

. . 9480 




®Time for the Sun to pass through 1 00 km of atmosphere at the constant 
rate Vj-qi (see text) . 
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TABLE 5.- DIVISION OF EARTH INTO ELEMENTS APPROXIMATELY EQUAL 
IN AREA TO A 5° X 5° BOX AT EQUATOR 


Latitude 

band, 

deg 

Area 

of 

band, 

kin^ 

AL,^ 

deg 

No.t> 

Area of 
each 
box, 
kni^ 

Cumulative 

area, 

km^ 

Percent of hemisphere 
area 

In band 

Cumulative^ 

0 5 

22 277 408 

5 

72 

309 408 

22 277 408 

8.72 

8.7 

5 10 

22 107 863 

5 

72 

307 053 

44 385 271 

8.65 

17.4 

10 15 

21 770 064 

5 

72 

302 362 

66 155 335 

8.52 

25.9 

15 20 

21 266 582 

5 

72 

295 369 

87 421 918 

8.32 

34.2 

20 25 

20 601 249 

5 

72 

286 128 

108 023 167 

8.06 

42.3 

25 30 

19 779 127 

6 

60 

329 652 

127 802 294 

7.74 

50.0 

30 35 

18 806 475 

6 

60 

313 441 

146 608 769 

7.36 

57.4 

35 40 

17 690 693 

6 

60 

294 844 

164 299 462 

6.92 

64.3 

40 45 

16 440 275 

6 

60 

274 004 

180 739 737 

6.43 

70.7 

45 50 

15 064 737 

8 

45 

334 771 

195 804 474 

5.89 

76.6 

50 55 

13 574 547 

8 

45 

301 656 

209 379 021 

5.31 

81 .9 

55 60 

11 981 046 

9 

40 

299 526 

221 360 066 

4.69 

86.6 

60 65 

10 296 362 

10 

36 

286 010 

231 656 428 

4.03 

90.6 

65 70 

8 533 317 

12 

30 

284 443 

240 189 745 

3.34 

94.0 

70 75 

6 705 328 

18 

20 

335 266 

246 895 073 

2.62 

96.6 

75 80 

4 826 307 

24 

15 

321 753 

251 721 380 

1 .89 

98.5 

80 85 

2 910 555 

40 

9 

323 395 

254 631 935 

1 .14 

99.6 

85 90 

972 653 

120 

3 

324 217 

255 604 588 

.38 

100.0 


^Longitude span of each box. 

*^umber of boxes in each latitude band, 
cprom equator to pole. 












TABLE 7.- MAXIMUM TANGENT LATITUDE AVAILABLE DURING SOLAR OCCULTATION 


MISSION AS A FUNCTION OF ORBIT INCLINATION AND ALTITUDE 


Altitude, 

km 

Maximum northern or southern tangent latitude, 
for specified orbit inclination, deg, of 

deg 

f 

a . 

^polar • 
deg 

0 

10 

20 

30 

40 

50 

60 

70 

80 

90 

200 

14.2 

24.2 

34.2 

44.2 

54.2 

64.2 

74.2 

84.2 

90 

90 

75.8 

400 

19.8 

29.8 

39.8 

49.8 

59.8 

69.8 

79.8 

89.9 

90 

90 

70.2 

600 

23.9 

33.9 

43.9 

53.9 

63.9 

73.9 

83.9 

90 

90 

90 

66.1 

800 

27.3 

37.3 

47.3 

57.3 

67.3 

77.3 

87.3 

90 

90 

90 

62.7 

1000 

30.2 

40.2 

50.2 

60.2 

70.2 

80.2 

90 

90 

90 

90 

59.8 

1200 

32.7 

42.7 

52.7 

62.7 

72.7 

82.7 

90 

90 

90 

90 

57.3 

1400 

34.9 

44.9 

54.9 

64.9 

74.9 

84.9 

90 

90 

90 

90 

55.1 

1600 

36.9 

46.9 

56.9 

66.9 

76.9 

86.9 

90 

90 

90 

90 

53.1 


^Minimum inclination for polar coverage at the specified altitude. 




TABLE 8.- NODAL PRECESSION RATES FOR CIRCULAR ORBITS AS A FUNCTION OF ALTITUDE AND INCLINATION 


deg/day, at altitude, h, km, of - 
200 300 400 500 600 700 800 900 1000 1100 1200 1300 1400 1500 


0 - 8.956 
10 - 8.820 
20 - 8.414 
30 - 7.752 
40 - 6.855 
50 - 5.749 
60 - 4.471 
70 - 3.057 
80 - 1.552 
90 0 

100 1.552 


- 8.495 
■ 8.366 
- 7.981 
• 7.353 
- 6.502 
- 5.454 
- 4.241 
- 2.900 
- 1 .472 
0 

1 .472 


- 8.064 

- 7.941 

- 7.576 

- 6.980 

- 6.172 

- 5.177 

- 4.026 

- 2.753 

- 1.397 

0 

1 .397 


- 7.661 - 7.283 
- 7.544 - 7.172 
- 7.197 - 6.842 
- 6.631 - 6.304 
- 5.864 - 5.575 
- 4.918 - 4.676 
- 3.825 - 3.636 
- 2.615 - 2.487 
- 1 .328 - 1 .262 
0 0 
1 .262 


- 6.929 - 6.597 - 6.285 - 5.992 - 5.716 - 5.456 

■ 6.824 - 6.496 - 6.189 - 5.900 - 5.629 - 5.373 

6.510 - 6.198 - 5.905 - 5.629 - 5.370 - 5.126 

- 5.710 - 5.440 - 5.187 - 4.948 - 4.723 

- 5.050 - 4.811 - 4.586 - 4.375 - 4.176 

- 4.236 - 4.035 - 3.847 - 3.670 - 3.503 

- 3.294 - 3.138 - 2.992 - 2.854 - 2.724 

- 2.252 - 2.146 - 2.046 - 1.952 - 1.863 

- 1 . 143 - 1 . 089 - 1.039 -.991 -.946 

0 0 0 0 0 
1.143 1.089 1.039 .991 .946 


1300 

1400 

1500 

- 5.211 

- 4.980 

- 4.762 

- 5.132 

- 4.904 

- 4.690 

- 4.896 

- 4.679 

- 4.474 

- 4.511 

- 4.311 

- 4.123 

- 3.989 

- 3.812 

- 3.646 

- 3.346 

- 3.198 

- 3.058 

- 2.602 

- 2.487 

- 2.378 

- 1.780 

0 

• 

1 

- 1 .626 

-.903 

-.863 

-.826 

0 

0 

0 

.903 

.863 

.826 












TABLE 9.- AVERAGE AND STANDARD DEVIATION OF g AS A FUNCTION OF 
LATITUDE BAND AND TIME FOR NOMINAL SOLAR OCCULTATION MISSION 


Latitude 

band, 

deg 

"q 

q 


t, 

days 

irst' 

days 

^last' 

days 

-90 -85 

68 

25.661 

4.136 

76.748 

74.498 

78.999 

-85 -80 

21 

26.182 

3.343 

73.759 

73.087 

74.430 


51 

27.936 

5.135 

80,746 

79.066 

82.426 

-80 -75 

135 

25.967 

4.250 

48.419 

43.917 

52,921 


15 

27.772 

6.288 

72.550 

72,080 

73,020 


59 

30.976 

5.559 

84.441 

82.493 

86.390 

-75 -70 

19^ 

34.62^ 

5.541 

241.426 

234,940 

247.911 


139 

31.52^ 

5.422 

308.910 

304,274 

313,546 


52 

28.515 

5.563 

42,136 

40.422 

43.850 


33 

28.656 

4.086 

54.063 

52,988 

55,138 


10 

29.391 

3.976 

71.711 

71.409 

72.013 


61 

30.486 

4.492 

66.473 

86.457 

90.489 

-70 -65 

62 

34.689 

4.439 

232.823 

230,772 

234,873 


7^ 

34.616 

6.072 

250,432 

247,978 

252,885 


A9 

36.057 

5.156 

302.594 

300.981 

304.207 


29 

34.241 

5.718 

314,554 

313.613 

315.494 


97 

34.029 

5.517 

348.097 

344,872 

351,323 


A1 

29.833 

5.235 

39.011 

37.667 

40,355 


20 

29.744 

4.644 

55.843 

55,205 

56,481 


7 

26.606 

4.740 

71,140 

70,939 

71.342 


58 

30.303 

5.520 

92.472 

90,556 

94,387 

-65 -60 

218 

35.534 

5.496 

205.142 

197,851 

212.434 


33 

35.927 

6,783 

229.630 

228.554 

230.705 


A7 

36.027 

5,837 

254.498 

252,952 

256.044 


39 

37.387 

6.080 

299.637 

298.361 

300,914 


20 

35.382 

4.862 

316.199 

315,561 

316.837 


41 

34.065 

5,096 

343,462 

342.118 

344.805 


56 

34.645 

5.166 

353.238 

351,390 

355,086 


36 

31.893 

4.433 

36,423 

35.247 

37,600 


13 

29.397 

4,648 

56,951 

56.548 

57.354 


5 

30.459 

4.002 

70,738 

70.604 

70.872 


50 

30.808 

5.477 

96.101 

94,454 

97,747 


60 

32.595 

6.032 

145.010 

143.028 

146,993 

-60 -55 

137 

35.189 

5,992 

193.214 

188.643 

197,784 


79 

36.715 

5,871 

215,123 

212.501 

217.744 


23 

37.276 

7,024 

227,747 

227.007 

228.487 


35 

37.409 

5,049 

257,254 

256,111 

258,397 


35 

39.774 

6.108 

297.151 

296.008 

298,293 


15 

38.535 

4,894 

317,374 

316,904 

317,844 


24 

37.441 

5.103 

341,278 

340,506 

342.051 


40 

37.003 

6,257 

356,464 

355.153 

357.774 


32 

33.717 

5,181 

34.138 

33,096 

35.180 


9 

33.551 

5.860 

57,689 

57,421 

57,958 


3 

— — -u 

32.648 

6,401 

70,470 

70.403 

70,537 




TABLE 9.- Continued 


Latitude 

band, 

deg 

"q 

q 


t, 

days 

^first' 

days 

^lastf 

days 

-60 -55 

41 

31,294 

4.090 

99,159 

97,814 

100.503 


123 

31,780 

5,306 

138.861 

134,760 

142.960 


153 

32,258 

5,034 

152,168 

147,060 

157,276 

-55 -50 

5<r 

34.305 

6.175 

186,795 

185,014 

188.576 


35 

38,288 

7,586 

218.954 

217.811 

220,097 


18 

36,945 

5,808 

226,368 

225,797 

226,940 


27 

40,473 

6.362 

259,338 

258.464 

260,212 


29 

40,560 

6.804 

295.000 

294.059 

295.941 


12 

37,695 

5,417 

316.261 

317,912 

318,650 


16 

35.772 

6.033 

339,868 

339,297 

3^0. ^39 


33 

34,741 

5,956 

358.917 

357.841 

359,992 


28 

31,785 

4.825 

32.122 

31,214 

33,029 


6 

26,395 

5.393 

58.193 

58.025 

58.360 


1 

27,824 


70,337 




32 

30,665 

5.205 

101,612 

100,570 

102.654 


42 

32,855 

5,174 

133.315 

131,937 

134,693 


58 

33,825 

5,797 

159,259 

157,343 

161,174 

-50 -45 

30 

35.435 

5,076 

183,972 

182,997 

184.946 


24 

34,538 

6.863 

220.938 

220.164 

221,711 


15 

36,578 

5,861 

225,259 

224,788 

225,730 


24 

39.186 

6,689 

261,052 

260.279 

261.825 


26 

38.398 

6.157 

293,152 

292.311 

293.992 


10 

36,522 

7,774 

319,019 

318,717 

319,321 


13 

37,143 

5.118 

338.827 

336.425 

339,230 


29 

32,028 

4.232 

361,000 

360.059 

361,941 


25 

30.718 

6,295 

30,340 

29.534 

31,147 


4 

32.588 

7,994 

56.528 

58,427 

58,628 


26 

31,345 

4,127 

103,561 

102,721 

104,402 


25 

32.399 

5,110 

131,063 

130,256 

131.870 


37 

33,629 

5,435 

162,452 

161,242 

163.662 

1 

1 

o 

19 

33,777 

5,801 

182.324 

181,719 

182,929 


33 

37.060 

7,693 

223,189 

221,778 

224,721 


19 

37,732 

6.154 

262,498 

261,892 

263.103 


23 

35,349 

5,176 

291,505 

290,765 

292,244 


7 

39,824 

5,275 

319,589 

319.338 

319.791 


10 

35.170 

7,414 

338,055 

337,753 

338,357 


25 

32.293 

4,718 

362,815 

362,008 

363,622 


23 

28,791 

4.686 

28,727 

27.988 

29.466 


2 

28,147 

4,969 

58,729 

58.695 

58,762 


20 

28.538 

4.848 

105,108 

104,469 

105.746 


18 

31,119 

5.477 

129,617 

129,045 

130.189 


28 

33,665 

5,775 

164,637 

163,729 

165,544 

-40 -35 

13 

33.440 

5,377 

181,248 

180.844 

181,652 


28 

34,141 

5.619 

223,456 

222,972 

223,997 


17 

34,652 

6.932 

263,708 

263,170 

264,246 


20 

34,019 

8.102 

290,059 

289.421 

290,698 


6 

31.961 

3.440 

320,025 : 

319.858 

320.193 


9 

34,182 

6,926 

337.418 

337.149 

337,686 


23 

32,007 

5,754 

364,428 

363,689 

365.168 




TABLE 9.- Continued 


Latitude 

band, 

deg 


-<f0 -35 



27.580 

35.396 

26.9A7 

29.780 
28.45<i 
32.161 
31.80<» 
32.202 

3^.016 

32.96^ 

36.362 

35.918 

32.121 

37.371 

29.274 

26.322 

27.517 

29.994 
28.232 
33.269 
32.098 
34.632 
32.029 
31.418 
33.317 
31.295 
28.423 
31.567 
27.217 
28.230 
26.433 
28.153 
27.936 
28.626 
23.828 
30.699 

30.780 
28.771 
28.483 

23.995 
27.004 
28.045 
26.577 
21.475 
29.946 
25.630 
30.739 
27.347 



t f 

days 

^first' 

days 

^last' 

days 

5.169 

27.251 

26.576 

27.920 


58.829 




70.339 



3«39^ 

106.352 

105.814 

106.690 

4.424 

128.541 

128.104 

128.978 

6.626 

166.351 

165.612 

167.091 

3.151 

172.047 

171.980 

172.114 

4.300 

180.508 

160.239 

180.777 

4.799 

222.534 

222.164 

222.904 

5.310 

224.502 

224*065 

224.939 

4.149 

264.784 

264.313 

265.255 

5.651 

288.782 

288.210 

269.353 

4.129 

320.394 

320.260 

320.526 

3.562 

336.915 

336.747 

337.082 


365.235 



4.628 

.699 

• 060 

1.338 

4.651 

25.903 

25.298 

26.508 

1.744 

70.441 

70.407 

70.474 

5.335 

107.360 

106.957 

107.764 

5.123 

127.700 

127.364 

128.036 

3.686 

167.763 

167.158 

168.368 

4.63<» 

171.274 

170.635 

171.912 

2.187 

180.003 

179.834 

180.171 

3.276 

221.761 

221.424 

222.097 

6.845 

225.376 

225.006 

225.746 

7.696 

265.726 

265.322 

266.129 

5.657 

287.639 

287.134 

268.143 

6.175 

320.696 

32 0.595 

320.796 

5.071 

336.512 

336.345 

336.680 

5.199 

2.043 

1.405 

2.682 

4.440 

24.659 

24.088 

25.231 

3.256 

70.610 

70.542 

70.677 

4.139 

106.166 

107.831 

106.504 

4.820 

127.061 

126.825 

127.296 

4.456 

169.502 

168.436 

170.568 

5.124 

179.632 

179.497 

179.767 

5.600 

221.021 

220.684 

221.357 

4.094 

226.184 

225.814 

226.554 

5.251 

266.566 

266.196 

266.936 

6.080 

266.630 

286.193 

287.067 

2.716 

320.930 

320.863 

320.997 

3.090 

336.177 

336.077 

336.278 

5.107 

3.321 

2.749 

3.692 

4.452 

23.516 

23.012 

24.021 

6.447 

58.792 

58.758 

58.826 

6.662 

70.812 

70.745 

70.879 

3.139 

106.841 

108.571 

109.110 

5.992 

126.556 

126.354 

126.758 

5.163 

169.502 

168.281 

170.722 


5 


S 








Latitude 

band, 

deg 


-20 -15 


-15 -10 



30.171 

29.709 
28.683 
25.135 

28.814 
28.184 
26.343 
25.990 
27.156 
25.586 
22.377 
29.820 
24.404 
25.969 
25.518 
27.529 
25.982 

12 I 26.756 

27.814 
26.847 
24.742 
29.093 
27.115 
26.358 
29.312 
27.961 
23.146 
28.917 
25.723 
24.913 

1 22.273 

12 23.663 

27.061 
24.942 
27.014 
26.511 
27.537 
25.347 

24.710 
23.646 
26.978 
26.172 
25.316 
25.336 
25.554 
23.986 
27.544 
20.448 
23.081 


7.279 

4.112 

5.066 

5.150 

4.036 

4.191 

1.803 

5.715 


1.908 

1.552 

3.470 

4.228 

3.452 

4.760 

.605 

3.103 

3.269 

4.709 

4.038 

I 

5.309 

5.084 

3.816 

5.604 

6.130 

3.439 

3.383 

4.911 

4.769 

4.573 

3.762 

3.634 

3.632 

4.416 

5.162 

3.177 

2.141 

5.837 


4.545 

5.624 

3.496 

3.517 

4.764 

5.744 


179.362 

220.247 

226.991 

267.306 

285.722 

321.097 

335.943 

4.464 

22.440 

58.657 

71.048 

109.379 

126.152 

167.743 

171.260 

179.193 

219.474 

227.798 

267.979 

284.661 

321.197 

335.776 

5.540 

21.398 

58.488 

71.351 

109.817 

125.848 

166.701 

171.966 

179.091 

218.700 

228.572 

266.618 

284.141 

335.676 

6.582 

20.389 

56.252 

71.668 

110.187 

125.612 

165.656 

217.859 

229.345 

269.190 

283.469 

335.610 


^first' 

days 


179.295 

219.877 

226.621 

267.004 

285.319 

321.064 

335.876 

3.960 

21.936 

58.623 

70.947 

109.177 

126.017 

167.272 

170.789 

179.160 

219.137 

227.428 

267.676 

284.511 

335.743 


20.927 

58.421 

71.217 

109.648 

125.747 

166.196 

171.798 

218.330 

228.235 

268.349 

283.839 

6.112 

19.918 

56.151 

71.554 

110.052 

125.545 

165.187 

217.456 

228.975 

268.955 

283.166 


^last» 

days 


219.070 

228.908 

268.888 

284.444 

7.053 

20.860 

58.353 

71.823 

110.322 

125.680 

166.129 

218.263 

229.715 

269.426 

283.771 
















TABLE 9.- Continued 


Latitude 

band, 

deg 

"q 

q 

-5 0 

14 

27.084 


4 

21.767 


6 

28.419 


3 

26.468 


3 

26,293 


15 

24,643 

0 5 

1 

26.336 


12 

24.839 


12 

26.808 


8 

24.739 


9 

23.967 


2 

27,405 


1 

26.141 


14 

24.431 


14 

24.846 


6 

27.674 


6 

26.466 


3 

21.474 


1 

19.184 


15 

24.406 

5 10 

2 

27.434 


14 

25.681 


12 

24,297 


7 

24,568 


8 

22.843 


2 

24.288 


1 

28.050 


14 

24.740 


14 

23,695 


5 

25.435 


8 

26.427 


3 

25,123 


1 

21.106 


16 

25,087 

10 15 

3 

24.352 


14 

26,452 


13 

25.588 


7 

25,015 


8 

25.517 


4 

23.850 


2 

2 7. ’9 60 


13 

24.719 


14 

25,164 


8 

31.358 


9 

25.104 


1 

24.555 


16 

27,646 

15 20 

3 

25.262 


16 

26,973 


3.287 
3.886 
A, 013 
<».128 
1.628 

3.833 

4.983 

i»672 

3.223 

4.288 

3.710 

4.812 

4.074 

3.433 

2.052 

4.870 

1.308 

4.288 

5.057 

6.132 

5.558 
.828 

3.515 

6.974 

4.640 

2.317 

4.441 

3.626 

2.969 

5.629 

2.558 
4.758 
4.235 
3.356 

5.554 

3.407 

4.060 


19. 

57. 

72. 
110 . 
125. 
164. 
179. 
217. 
230. 
269. 
282. 
321. 
335. 

8 . 
18.472 
57.645 
72.463 
110.659 
125.274 
163.641 
179.190 
216.144 
230.960 
270.234 
282.257 
321.025 
335.679 
9 
17 


^first' ^last' 
days days 


977 19,851 

881 58.083 

890 72,227 

389 110.524 
342 125,477 
165.120 

217.388 
230.523 
269.964 
283,099 
321,195 


9 


18.910 

57.814 

72.631 

110.727 


9,9 


17,968 

57,410 

73,169 

110.929 

















TABLE 9.- Continued 


Latitude 

band, 

deg 


15 20 




q 

28 

• 884 

25 

.955 

23 

• 981 

24 

• 528 

25 

.460 

27 

.616 

26 

.950 

28 

.301 

30 

.435 

32 

.529 

31 

.165 

28 

• 910 

27 

.365 

26 

.412 

27 

.431 

24 

.541 

26 

.027 

27 

.957 

21 

.097 

28 

.919 

27 

• 692 

28 

.550 

30 

• 849 

27 

.251 

29 

.245 

30 

• 381 

26 

.703 

30 

• 008 

26 , 

,579 

24 

.528 

26 , 

,907 

27 , 

,893 

28 

,827 

30 , 

,484 

29 , 

,347 

34 , 

,970 

30 . 

,896 

29 , 

,499 

33 , 

,529 

28 , 

,274 



26.381 

28.961 

26.401 

28.296 

29.480 

31.005 

33.124 

38.229 


3.100 

5,434 

2.322 

3,140 

2,009 

4.544 

3.076 

3.886 

4.675 


4,013 

1.420 

4.238 

4,823 

3,514 

3.955 

4,098 

1.283 

4,647 

4,389 

6,598 

4.346 

.643 

4.640 

2,764 

4.683 

2.620 

1.696 

4.988 

4,410 

4,405 

4.363 

3.443 

5.271 

4,206 

4.701 

3,545 

4,572 

5.326 

1,784 

6,926 

5,278 

5.077 

5.069 

6,306 

5,142 

2.012 


t, 

days 


232,709 

271.176 

281,214 

320,553 

335.983 

11.324 

15,648 

56,366 

74,145 

111.065 
125.205 
160.413 
179.826 
213.084 
233,684 
271.614 
280.743 
320,216 
336.287 

12.299 
14.639 
55,794 
74,885 
125,305 
159.236 
180.194 
211 ,908 
234,726 
272,032 
280.272 
319,711 
336.657 
13.425 
55,155 
75,726 

111.066 
125,472 
157.959 
180.630 
210,597 
235,869 
272,389 
279,834 
319,072 
337,162 

13,339 

54.381 

76,701 

110,967 


tfi 

rst' 

^last' 

days 

days 

232^ 

,271 

233 , 

.146 

270 i 

» 97 ^ 

271 , 

.378 

281 , 

.012 

281 , 

,416 

320 , 

,452 

320 , 

,654 

335 < 

.882 

336, 

,084 

10 , 

,887 

11 , 

,761 

15 , 

,210 

16 , 

,085 

56 , 

,130 

56 , 

,602 

73 , 

1042 

74 , 

,448 

159 , 

875 

160 , 

951 

179 , 

,692 

179 , 

^960 

212 . 

546 

213 , 

,622 

233 . 

213 

234 . 

155 

271 , 

446 

271 , 

782 

280 . 

541 

L 80 . 

945 

320 , 

048 

320 , 

385 

336 « 

152 

336 . 

421 

11 . 

828 

12 . 

770 

14 , 

134 

15,143 

55 , 

525 

56 . 

063 

74 , 

515 

75 , 

255 

125 , 

271 

125 , 

338 

158 . 

665 

159 . 

808 

180 , 

027 

160 . 

362 

211 . 

336 

212 . 

4 79 

234 . 

222 

235 . 

231 

271 , 

850 

272 . 

186 

280 . 

070 

280 . 

474 

319 , 

442 

319 , 

980 

336 , 

489 

336 . 

825 

12 . 

837 

14 . 

067 

54 , 

852 

55 . 

457 

75 . 

322 

76 , 
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405 

125 . 
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157 . 
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158 . 
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429 
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831 

209 . 
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211 . 

269 

235 . 
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272 , 
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272 . 

523 

279 , 

666 

280 . 

002 

318 . 

769 

319 . 

375 

336 . 

893 

337 , 

431 

11 . 

781 

14 . 

922 

53 , 

977 

54 , 

785 

76 , 

197 

77 , 

205 

110 . 

933 

111 . 

000 














TABLE 9.- Continued 


Latitude 

band, 

deg 

"q 

q 


t, 

days 

^first' 

days 

^last' 

days 

30 35 

5 

30.783 

4.531 

125.740 

125.606 

125.874 


21 

31.931 

5.971 

156.580 

155.908 

157.253 

35 

e 

30.199 

4.152 

181.133 

180.898 

181.368 


22 

30.484 

5.241 

209.151 

208.445 

209.857 


19 

27.340 

4.389 

237.113 

236.508 

237.718 


6 

30.168 

5.303 

272.759 

272.591 

272.927 


6 

30.535 

2.442 

279.430 

279.262 

279.598 


15 

27.836 

5.840 

318.231 

317.760 

318.702 


12 

27.022 

5.108 

337.869 

337.498 

338.239 


22 

31.420 

5.511 

11.007 

10.301 

11.713 


18 

31.286 

5.473 

15.561 

14.989 

16.132 


15 

32.777 

4.649 

53.439 

52.969 

53.910 


18 

36.612 

5.379 

77.844 

77.273 

78.416 


3 

32.479 

1.109 

110.800 

110.733 

110.867 


5 

32.031 

6.972 

126.075 

125.941 

126.209 


22 

32.807 

4.917 

155.135 

154.429 

155.841 

40 45 

11 

32.341 

5.991 

181.771 

181.435 

182.106 


25 

31.518 

5.015 

207.572 

206.765 

208.378 


21 

27.280 

5.830 

238.456 

237.786 

239.130 


5 

25.687 

5.318 

273.129 

272.995 

273.264 


6 

27.315 

7.262 

279.026 

278.858 

279.194 


19 

27.810 

3.835 

317.088 

316.483 

317.693 


18 

29.438 

4.143 

338.878 

338.306 

339.449 


28 

35.174 

5.831 

9.326 

8.419 

10.234 


21 

33.237 

6.271 

16.872 

16.200 

17.545 


18 

36.353 

5.461 

52.330 

51.758 

52.901 


21 

36.644 

5.817 

79.155 

78.483 

79.827 


4 

32.576 

6.438 

110.565 

110.465 

110.666 


8 

35.013 

3.800 

126.51l| 

126.276 

126.745 


26 

32.968 

5.566! 

153.522 

152.661 

154.362 

45 50 

13 

34.420 

5.722 

182.576 

182.173 

182.979 


29 

32.292 

4.413 

205.757 

204.816 

206.698 


24 

28.522 

5.728 

239.971 

239.197 

240.744 


6 

26.848 

2.706 

273.500 

273.332 

273.668 


5 

27.186 

3.082 

278.655 

278.521 

278.790 


27 

29.294 

4.930 

315.541 

314.667 

316.415 


28 

28.766 

5.357 

340.424 

339.516 

341.332 


39 

33.522 

4.396 

7.074 

5.797 

8.352 


26 

33.004 

7.134 

18.452 

17.612 

19.293 


24 

37.124 

5.732 

50.918 

50.144 

51.691 


25 

38.394 

6.348 

80.701 

79.895 

81.508 


6 

38.470 

5.673 

110.230 

110.063 

110.398 

1 

10 

39.183 

9.846 

127.115 

126.813 

127.417 


27 

3 7 .3 83 

6.616 

151.740 

150.866 

152.614 

50 55 j 

17 

34.507 

6.490 

183.583 

183.046 

184.121 


32 

34.118 

4.932 

203.707 

202.665 

204.748 


27 

30.671 

5.148 

241.685 

240.811 

242.559 


5 

30.377 

5.857 

273.870 

273.736 

274.005 


6 

32.754 

6.227 

278.285 

278.117 

278.453 









TABLE 9.- Continued 


Latitude 

band, 

deg 

"q 

q 


days 

^first' 

days 

^last' 

days 

50 55 

AO 

29.064 

5.337 

313.289 

311.979 

31A*600 


55 

32.630 

5.581 

343.214 

341.399 

345.029 


68 

34.518 

5.087 

3.478 

1.226 

5,730 


35 

35.413 

5.828 

20,503 

19.360 

21.646 


31 

40.290 

7.053 

49,069 

48,060 

50,077 


30 

38.366 

5.370 

82,550 

81,575 

83,524 


7 

41.417 

6.598 

109,794 

109.593 

109.996 


13 

38.759 

5#616 

127,887 

127.484 

120.209 


32 

35.471 

6.216 

149,757 

148,716 

150.799 

55 60 

22 

34.393 

4.199 

184.893 

184.186 

185.598 


40 

34.246 

6.051 

201.302 

200.044 

202,598 


31 

32.032 

5.789 

2A3#63A 

2A2*626 

2AA.6A2 


6 

30.978. 

1.494 

274,241 

274.073 

274.409 


5 

26.081 

3.566 

277,915 

277.780 

278,049 


59 

30.535 

5.703 

309.962 

306.013 

311.911 


300 

33.357 

5.643 

355,144 

3A5.096 

365.191 


18 

34.886 

6.805 

.588 

.016 

1.159 


55 

39.732 

6.376 

23,526 

21,713 

25.343 


49 

40.755 

6.403 

46.380 

44.767 

47,993 


3A 

41.275 

8.025 

84,701 

63.592 

85,810 


9 

42.032 

6.570 

109,258 

106.989 

109,526 


18 

37.177 

6.765 

128,927 

128,357 

129.498 


36 

36.468 

6.571 

147,472 

146,296 

148.648 

60 65 

35 

37.354 

5.640 

186,807 

185,665 

187,949 


51 

33.732 

6.151 

198,297 

196,617 

199,977 


35 

30.699 

4.546 

245.852 

244,710 

246.995 


5 

32.995 

4.785 

274,611 

274.477 

274,746 


5 

33.019 

2.253 

277.578 

277,443 

277,712 

. 

ee 

29.616 

4,651 

305,023 

302.099 

307.946 


280 

39.655 

6.887 

35.055 

25.410 

44,699 


39 

40.385 

8,014 

87,154 

85,877 

88,431 


12 

43.536 

7.694 

108.553 

106.184 

106.922 


27 

40.856 

6.654 

130.436 

129,565 

131.312 


AA 

39.316 

7.237 

144,784 

143,339 

146.229 

65 70 

128 

37.071 

6.395 

192,283 

168.016 

196.549 


A2 

30.776 

5,562 

248.440 

247,062 

249,818 


7 

29.283 

1.888 

275.015 

274.813 

275,217 


6 

28.247 

4,050 

277,207 

277.039 

277,376 


113 

28.685 

4,556 

298.269 

294.506 

302,032 



43.460 

7,941 

89.943 

68.496 

91.387 



40.793 

5,012 

107,613 

107,109 

108.116 



41.352 

6,510 

133,193 

131,379 

135,007 


71 

41.087 

7,592 

140.920 

138.568 

143,272 

70 75 

48 

32.906 

6,217 

251,411 

249.885 

252,976 


16 

29.082 

4.904 

276.043 

275.285 

276,971 


mnm 

27.617 

4,136 

290,205 

285,971 

294,438 



42.899 

7,144 

93,067 

91,455 

94.680 



43.501 

7,304 

106,337 

105,632 

107.042 


Bil 

43.508 

7,943 

136,787 

135,074 

138,501 











Latitude 
band , 
deg 


75 80 



56 33,257 

155 26,908 

59 <» 5,621 

36 <» 5 , 85 ^ 

67 33,024 

9 26,304 

67 47,808 

213 26,542 


5,522 

4,465 

8,094 

6,804 

5,819 

5,352 

7,041 

5,081 


254.892 

280,756 

96,696 

104,389 

259.025 

275,894 

100,929 

268,433 


^first' 

^last» 

days 

days 


253,044 

275.958 

94,747 

103.213 

256,807 

275,625 

98.712 

261.310 












TABLE 10.- LEAST SQUARES PARAMETERS AND RELATED QUANTITIES FOR MODELING 
BEHAVIOR OP q WITHIN A LATITUDE BAND 


Latitude 


Model parameters 



band 

, deg 

"q 




^ 2 

% 

n a 

a 






Q. 

A 

y, deg 



-90 

-85 

1 






-85 

-80 

2 






-80 

-75 

3 

63.84 

38.03 

-142.9 


100.0 

-75 

-70 

6 

31 .86 

3.41 

-115.9 

0.05 

99.4 

-70 

-65 

9 

32.07 

3.54 

175.8 

1.70 

87.4 

-65 

-60 

12 

33.67 

3.24 

-172.8 

.61 

92.5 

-60 

-55 

14 

35.30 

3.75 

163.1 

.28 

96.8 

-55 

-50 

13 

34.35 

5.60 

-172.4 

2.46 

89.7 

-50 

-45 

13 

34.38 

3.49 

-172.6 

2.03 

78.0 

-45 

-40 

12 

33.34 

5.02 

-168.5 

2.72 

84.6 

-40 

-35 

14 

32.02 

2.50 

-168.9 

5.19 

40.9 

-35 

-30 

15 

32.09 

3.45 

-175.7 

5.11 

54.2 

-30 

-25 

13 

29.73 

2.88 

-146.8 

3.03 

62.2 

-25 

-20 

14 

27.40 

1 .40 

-127.6 

8.37 

13.3 

-20 

-15 

15 

26.89 

1.45 

-149.0 

4.66 

21 .3 

-15 

-10 

15 

26.85 

.62 

32.1 

3.27 

7.1 

“10 

-5 

13 

25.53 

.88 

115.6 

2.60 

15.3 

-5 

0 

12 

25.06 

1 .25 

-59.6 

6.06 

13.0 

0 

5 

14 

24.90 

1.35 

138.2 

5.16 

18.3 

5 

10 

14 

24.91 

.08 

85.4 

3.69 

.1 

10 

15 

13 

25.98 

.69 

12.4 

4.63 

6.3 

15 

20 

14 

27.65 

2.85 

-18.9 

2.54 

67.2 

20 

25 

13 

27.34 

2.11 

-11.8 

4.29 

38.3 

25 

30 

13 

28.99 

2.69 

-16.8 

3.03 

62.8 

30 

35 

13 

30.20 

3.66 

-17.0 

4.56 

67.6 

35 

40 

14 

30.92 

2.66 

-4.5 

2.96 

60.4 

40 

45 

14 

I 31.64 

4.52 

6.5 

2.38 

84.6 

45 

50 

14 

33.16 

5.95 

-11.8 

.49 

97.9 

50 

55 

14 

34.93 

4.68 

-1.2 

2.42 

85.4 

55 

60 

14 

35.32 

5.99 

2.4 

3.09 

88.2 

60 

65 

11 

36.52 

5.29 

-2.2 

2.74 

89.7 

65 

70 

9 

35.49 

6.76 

-17.6 

.96 

98.2 

70 

75 

6 

35.83 

8.45 

-41 .0 

.66 

99.3 

75 

80 

4 

36.19 

14.90 

-59.0 

.89 

99.7 

80 

85 

3 

37.10 

17.57 

-61.8 


100.0 

85 

90 

1 







^Arbitrary but consistent units. 
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Figure 1 Orbit geometry for viewing sunrises and sunsets from a spacecraft, 


A 

r 



N = (sin sin i, -cos sin i, cos i) 

H = (N X r) 

Figure 2.- Definition of pointing angles to the Sun relative to a spacecraft. 






60 



Time from launch, days 


(a) Earth tangent latitude versus time from launch (see fig. 1). 



Time from launch, days 


(b) Pointing angle P versus time from launch (see fig. 2) . 

Figure 4.- Orbital data for a 1-year solar occultation mission 
using a 70°, 600-km circular orbit. 



deg 


deg 



0 30 60 90 120 150 180 210 240 270 300 330 360 

Lp deg 

(c) Earth tangent latitude versus longitude for sunrises. 



0 30 60 90 120 150 180 210 240 270 300 330 360 

Ly. deg 

(d) Earth tangent latitude versus longitude for sunsets. 

Figure 4.- Continued. 
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Time from launch, days 

(e) Apparent vertical velocity of the Sun's image relative to the horizon 
at the instant of sunrise or sunset. t-|oo “ Time to pass through 
100 km at constant rate Vj-gi. 


Figure 4.- Concluded. 



-24 



0 2 4 6 8 10 12 14 

Time from 100 km, min 


Figure 5.- Sun viewing requirements and surface coordinates during a solar sun 
set and sunrise near the beginning of a period of occultation (day 7 of the 
nominal mission). h 0 = Altitude of center of Sun's image above the local 
horizon. 


+ Sunset 


• Sunrise 



(b) During first 30 days after start of occultations. 

Figure 6.- Tangent latitudes and longitudes of sunrises 
and sunsets during nominal mission. 




+ Sunset . Sunrise 



(c) During first 90 days after start of occultations. 
Figure 6.- Concluded. 
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Figure 7.- Distribution of occultation measurements as a function of latitude and longitude boxes of 
approximately equal area (3 x 10^ km^) for the nominal solar occultation mission. 




(a) Viewing periods of continuous sunlight given a 
value of solar declination. 



(b) Continuous sunlight never occurs. 



(c) Continuous sunlight is possible regardless of solar declination. 

Figure 8.— Examples of limiting cases and geometries for viewing 
sunrises and sunsets from spacecraft orbits. 
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Figure 9.- Orbit and solar declination conditions required to produce 
periods of continuous sunlight during solar occultation missions. 



Figure 10.- Geometry for establishing the tangent latitude 
of a solar occultation measurement. 




n - RAg, deg 


(a) i = 50 °. 

during solar occultation missions for 50 ° and 70 °, 
plane orientation relative to the Sun, with solar 





o 



- RAg, deg 


(b) i = 70°. 


Figure 11.- Concluded 
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Figure 13.- Envelopes of tangent latitudes available during solar occultation 
missions at 600 km and various inclinations. 
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Figure 15.- Distribution of measurement (tangent) latitudes for 1-year simulated 
solar occultation missions in 600-km circular orbits. 
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